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PREFACE 


The ob jective of this work was to perform a study, prepare a Minitrack error 
model with as many of the coefficients as practical being established by pre-flight 
calibration measurements, prepare a computational method to utilize the error model,, 
and to monitor the MOTS camera drive stability using a diffraction grating. 

The NAP-II program, which has the capability of simultaneously solving for 
orbits and tracking station error model terms, was modified to substantially increase 
its computational speed and substantially reduce its disk memory requirements, thus 
making it possible to meet the contract objective within the required time-scale. 

Twenty-five short optically determined reference arcs were used to calibrate 
the Fort Myers Minitrack station. 

The NAP program was used to obtain a simultaneous four-arc solution (total 
length 25 days) of the 4 orbits and Minitrack station calibration numbers (for all 
stations) . 

A diffraction grating was designed and used for monitoring the stability of the 
Fort Myers MOTS camera drive . 

The multi-arc approach to the self-calibration of Minitrack stations appears 
to give very good results when judged on the comparison between Minitrack and 
optically determined orbits. The four arcs processed under this contract showed RMS 
position differences of 67m, 86m, 124m and 168m, respectively. This compares 
favorably with an RMS position difference of 165m for the first arc based on "aircraft 
calibrated" Minitrack measurements. 

Use of the diffraction grating apparatus developed under this contract provides 
a simple and inexpensive means of correcting plate errors caused by the (periodic) 
instability of thu MOTS camera drive. 

It is recommended that the approach to Minitrack self-calibration developed 
under this contract be further tested using other satellites (than GEOS-I) before 
"aircraft calibration" is finally abandoned. 


in 
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SECTION 1 


MONITORING OF MOTS STABILITY 
BY MEANS OF A 

DIFFRACTION GRATING APPARATUS 
1.1 INTRODUCTION 

The NASA MOTS camera and the U. S. Air Force PC- 1000 camera employ the 
same 1000mm f/5 Baker telephoto lens. The PC-1000 uses an alta-azimuth mount and 
is locked in a fixed orientation (relative to the earth) throughout each exposure. A 
precisely timed shutter provides intermittent exposures of star trails to serve as control 
points. The MOTS camera, by contrast, uses an equatorial mount and is sidereally 
driven to maintain a fixed orientation with respect to the right ascension declination 
frame. By virtue of this mode of operation, MOTS does not require a precisely timed 
shutter and provides a greater abundance of stellar images than a PC- 1000. Both MOTS 
and PC-1000 can potentially produce accuracies of about 0.6 seconds of arc for satellite 
directions. This potential, however, may not be routinely realized. 

One of the key advantages claimed for the PC- 1000 (and for the fixed camera 
mode of operation, in general) is that any significant drift in the orientation of the 
camera throughout an exposure can be detected (and generally corrected for) by virtue 
of separate reductions performed on each sequence of stellar exposures. Brown, 
(Reference 1) reports results of an investigation of the stability of a PC-1000 over a 
period of about one half hour. Reductions were performed on exposures made at 5 
minute intervals, leading to results plotted in Figure 1, which depicts the temporal 
variation of the three angular elements of orientation about their respective means. 

The plotted results for hour angle are normalized by the customary process of multipli- 
cation by cosine of declination. Each point is accompanied by a vertical bar defining 
its plus and minus 1 sigma confidence intervals. The 1 sigma values for declination and 
normalized hour angle are slightly less than 0.2 second of arc; for swing angle, they are 


1 




generally about ten times greater, averaging close to 2 seconds of arc. This disparity in 
sigmas is attributable to the fact that the focal length (1000mm) of the PC- 1000 is about 
10 times greater than the semi-diagonal of the plate format. The projective effect of 
an error of 0.2 second of arc is the direction of the camera axis is equivalent to that of 
an error of aboui 1 micron on the plate. By the same token, the projective effect of an 
error of 2 seconds in swing angle is equivalent to that of an error of about 1 micron near 
the edge of the plate. Thus, there is no actual projective disparity in the relative sigmas 
of the angular elements. 

Figure 1 clearly demonstrates that significant changes in the orientation of the 
PC- 1000 can occur throughout a period as short as 5 minutes. Indeed, changes in hour 
angle and declination of as much as 2 seconds of arc in 5 minutes are not unusual. Be 
this as it may, the fact is that in a vixed camera operation, any changes of significance 
can be detected and their effects on directions of satellites can be largely removed by an 
interpolative process. Thus while stability is a problem in PC-1000 operations, it is a 
problem that can be routinely overcome. 

In general, MOTS is subject to the same sources of instability as the PC-1000 
(e.g., thermal imbalances, wind loading, etc.). However, MOTS is also subject to 
drifts induced by the sidereal drive. Whatever their origin, small instabilities in MOTS 
orientation cannot be routinely detected because their primary effect is merely to cause 
a slight enlargement in the sizes of stellar images. This is a matter of no consequence 
insofar as the plate reduction itself is concerned, for instability affects all stars alike 
and hence does not degrade the residuals produced by the least squares adjustment. 
However, flashes are affected by the instantaneous rather than the integrated effects of 
instability. Hence, satellite directions derived from MOTS observations inherit fully ' 
the instantaneous departures in orientation of the camera from the mean orientation 
deduced from the stellar exposures. Clearly then, a simple, practical, routine means 
for monitoring the stability of MOTS orientation tl.roughout an exposure would constitute 
a significant advance. The objective of the present investigation is to evaluate an 
approach that holds promise in this regard. 


3 



1 .2 GENERAL CHARACTERISTICS OF ERROR INDUCED BY MOTS SIDEREAL DRIVE 


A study reported by Harris, Cartwright and Oosterhout (Reference 2) provides a 
good understanding of the nature of the errors induced by the MOTS sidereal drive. The 
data analyzed in Reference 2 consists of MOTS images of a stationary collimator generating 
an artificial star produced by a pinhole illuminated by a strobe light of 1 millisecond dura- 
tion flashing at 10 second intervals. The exposures of the flashing collimator were made 
over a period of 16 minutes with the MOTS axis in a nominally horizontal orientation and 
with the camera being driven in the normal sidereal mode. This generated a set of almost 
100 successive images which, with a flawless drive, would ideally be at a constant declina- 
tion and spaced at equal increments of hour angle. Departures of actual images from their 
ideal positions provide the desired measure of errors of the drive. In the analysis provided 
in Reference 1, polynomials in time were fitted to the measured x and y coordinates of 
the collimator images in order to account for low order effects such as the slight curvature 
of the trials. The residuals about the fitted polynomials reflect the combined effect of 
plate measuring errors and higher order sidereal drifts (slow thermal drifts and the like would 
be absorved by the fitted polynomials). 

Graphs of residuals in right ascension and declination for a particular trial presented 
in Reference 2 are reproduced in Figures 2 and 3 below. The most significant finging of 
Reference 1, namely the existence of a well defined periodic drift in right ascension, is 
well illustrated in Figure 2. The period of the drift is almost precisely 90 seconds and is 
attributed to the interaction of the worm and sector gears of the drive system. The ampli- 
tude of the drift is reported in most cases to range between 2 and 3 arc seconds in right 
ascension and to be only about half as great in declination. 

Because of geometrical considerations, the procedure employed in Reference 1 to 
monitor MOTS stability is not one that can be practically adopted to effect a routine 
calibration of drift affecting operational plates. The procedure has served well, however, 
to define the essentia! characteristics of MOTS drift, thereby making clear the problem to 
be overcome by any method designed for routine operational calibration of drift. 
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SECONDS 


FIGURE 3 


-Sudbury Data Pass No. 6 Quadratic Residuals A Y vs. Time 


. Angular stability of MOTS in declination as determined by flashing collimator 
method reported in reference [2] (5 microns = 1 second of arc). 







1.3 THE DIFFRACTION GRATING METHOD 

The present Investigation is concerned with the evaluation of a method of monitoring 
MOTS stability proposed to NASA by DBA Systems. As originally conceived, it involved 
placement of a coarse diffraction grating in front of the MOTS lens to generate measurable 
first order diffraction images of sufficiently bright stars. A series of six exposures would be 
made at suitable intervals with orientation of the grating being changed in nominal 15 
degree steps between exposures. This sequence of exposures was to provide material for a 
'precalibration 1 of MOTS drift. Shortly before the exposure of the satellite, the diffraction 
grating would be automatically retracted; immediately following the completion of the 
satellite observations, the grating would be returned in front of the objective for a series 
of six 1 postcalibration 1 exposures at 15 degree steps. A ten times enlargement of the image 
generated by this process for a sufficiently bright star would have the general appearance 
shown here: < O j . Now, in theory, the mean of a pair of opposing diffraction images 

V / 

will coincide precisely with the center of the associated primary central image. It follows 
that, should the central image wander slightly in position (because of drift of the mount), 
the locus of the wander can be reconstructed from the means of diffraction images taken at 
different times. The principle of the diffraction grating approach is thus straightforward 
and simple. 

From early experiments, a more satisfactory operational method of implementing 
the diffraction grating approach evolved. It was learned that it was not really necessary 
to employ the entire aperture of the MOTS in order to obtain satisfactory images; rather, 
a diffraction grating having a diameter of less than half that of the MOTS objective could 
provide altogether acceptable results. This finding meant that the grating could remain in 
place and undergo its series of step rotations throughout the period in which satellite obser- 
vations were being made. This obviates the need for interpolation between pre and post- 
calibration series of exposures and provides instead a running calibration spanning the 
interval of actual interest. 
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The final version of the apparafus that - evolved from exploratory exercises is shown 
at half scale in Figure 4. The grating itself is of 3. 5- inch diameter and is supported over 
the center of the field of the objective by a three-legged spider. The grating and spider 
obscure about 20 percent of the area of the aperture. With the particular grating employed, 
25 percent of the light incident on the grating goes into the central image. Accordingly, 
the effective loss of light due to the imposition of the grating is only about 15 percent. 

The grating is rotated by means of a belt-driven turntable connected to a motor. A remote 
manually operated, pushbutton switch causes a spring loaded solenoid plunger to retract 
from one of a series of slots spaced at 15 degree intervals near the rim of the turntable. 

This retraction immediately trips a microswitch which starts the drive motor. The turntable 
is driven 15 degrees, whereupon the spring loaded plunger slips into the next slot, tripping 
a microswitch which turns off the drive motor. This rotation of 15 degrees is accomplished 
in about one half second. Thus the operator merely activates the pushbutton switch for an 
instant whenever he wishes to effect an increment of rotation of the grating* 

The grating itself is contructed from nylon monofilament of 0.015 inch diameter 
spaced at 0.030 inch intervals (thus the width of the grating openings is equal to the 
diameter of the monofilament, causing half the light to be absorbed by the grating). This 
generates about N = 115 openings over the 3.5-inch aperture of the grating. A Xerox 
print of the grating is provided by Figure 5. For the adopted design, the theory of diffrac- 
tion gratings shows that one fourth of the light impinging on the grating will be directed to 
the central image; each of the first order images will receive about 10 percent of the 
original light (or, more precisely, a fraction of 1 /V s of the light). The angular distance 
between the central image and each of the first order images is 8 = x/d where X is the 
wavelength of the light and d is the center to center spacing of the grating divisions. For 
blue light (X = 0.0005mm) and a spacing of d = 0.030 inches (or 0.75mm), the value of 0 
becomes 1/1500 radian. Because the MOTS focal length is nominally 1000mm, this corre- 
sponds to a spacing of about 0.670mm on the MO T S plate. 
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FIGURE 4. Half scale drawing of diffraction grating apparatus employed to monitor 
MOTS stability. 9 
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1.4 EXPERIMENTAL PROCEDURE 

Short visits were made to the MOTS camera at the Ft. M/ers STADAN station in 
April, September and November 1970 to gather data needed for the evolving design of 
the apparatus just described. Finally, on the night of December 2, a highly successful 
series of observations was made with the final version of the diffraction grating. In this 
section and the next we shall review the results of this test. 

Because of other tests being conducted at the same time, the MOTS cameras had 
been refocussed for optimum imagery with blue light. For this reason, an orthochromatic 
emulsion (103 Ao) was employed for the diffraction grating experiment. For the data 
plate that was actually reduced, the MOTS was pointed toward the star Persius, the 
drive was started at 01 h 55“ EST (December 3), and the initial exposure was started one 
minute later. Thereafter, exposures were made in accordance with the pre-established 
log reproduced on the next page. The total exposure for a complete set of 12 successive 
grating positions amounted to 180 seconds (or to two complete periods of the drive error), 
the grating being shifted at 15 second intervals throughout the exposure. Thus each 
central image of a star received the full 180 second exposure (not only from the grating 
but also from the unobscured aperture of the objective, as well), while each pair of 
diffraction images received a 15 second exposure (from the grating only). To simulate 
the taking of multiple plates, this process was repeated a total of six times, the drive 
being turned off for one minute between each sequence. Although all six sequences 
were recorded on a common plate, the practical end result was equivalent to what would 
have been obtained from the exposure of a single sequence on each of six separate plates. 
For this reason, we shall hereafter refer to the results from exposure sequences, 1 through 
6 as being from plates 1 through 6. 

The star £ Persius toward which the camera was initially aimed, generated diffrac- 
tion images that were too large for precise measurement. However, usable diffraction 
images were generated by eight other stars ranging in stellar magnitude from 2.5 to 4.0. 
Inasmuch as a set of diffraction images from a single star is sufficient for monitoring 
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stability, there was ample data from which to make a selection for reduction. In order 
to obtain a measure of redundancy, we selected two stars from each 'plate' for measure- 
ment and reduction. A 40X enlargement of the diffraction images generated by one of 
the two selected stars is provided in Figure 6. The total exposure of the central image 
is almost 600 times greater than the exposure of a given diffraction image. This explains 
the relatively large diameter of the centrl image. Although somewhat ragged because of 
their small size (about 20 to 30 microns on the plate), the diffraction images are, none- 
theless, of a satisfactory quality for precise measurement. 

1.5 EXPERIMENTAL RESULTS 

On each of the 6 'plates', the 24 diffraction images (12 pairs) for each of the two 
selected stars were measured by an experienced operator on DBA's Mann comparator. 

The first phase of the data reduction consisted simply of determining the mean of the 
coordinates of each of the 12 pairs of diffraction images and then of subtracting from 
these the grand mean obtained by averaging the set of 12 means. If there were no drift 
of the camera and no measuring error, the departures from the grand mean would be 
precisely zero for each averaged pair of diffraction images. Significant, systematic 
departures were found to exist, and, within acceptable measuring tolerances, were 
found to be the same for the two stars measured on each plate. Accordingly, the depar- 
tures for the two stars measured on each plate were themselves averaged for corresponding 
pairs of diffraction images. 

The first stage in the analysis of the results consisted merely of plotting as functions 
of time the x (right ascension) and y (declination) departures obtained by the process just 
described. In addition, a program was written to perform a least squares regression based 
on functions of the form: 

' _ ■ , ' . 2tt 2tt 

6x 3 -ap + r 3 + a 2 sin r 3 + a 3 cos r 3 


6y i - b 0 + b x T j + bg sin 


2tt 

90 f i 


+ b 3 cos 


2t t 

90 T J 
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in which. 


6xj , 6y 3 - 


observed departures in x and y from grand mean of averaged 
coordinates for i th pair of diffraction images. 


Tj 


ao,b 0 

a i/bi 

a 2/ a 3 

^ 2/^3 


= relative time of midpoint of exposure of ! th pair of diffraction 
images 

= 15 (j — 1 ) second (i.e., r i = 0/ r 2 = 15,. . . ,Tis = 155 seconds), 
= zero set coefficients,, 

= linear drift coefficients, 

= periodic drift coefficients (90 second period assumed). 


Values of the observed 6Xj , 6 y j plotted against r i are presented in Figures 7 through 12. 
Also plotted in the same figures are corresponding values of 6x 3 , 6 y j as computed from 
the fitted functions. 

In reviewing Figures 7 through 12, we find that the departures in right ascension 
and declination are not randomly distributed about zero, but rather, for the most part, 
display fairly well defined trends. Flowever, with a few notable exceptions, the fit of 
the regression functions can be said to be only fair. This is especially clear from the 
results presented in Table 1 on the following page. These show that goodness of fit 
indicated by S x and S y is generally only modestly better than the dispersion about zero 
means as indicated by s x and s y . Indeed, in a few instances, values of S x and S y are 
slightly larger than their counterparts s x and s y , a result attributable to the fewer degrees 
of freedom associated with the residuals from the fitted functions (i.e. , 8 degrees of free- 
dom are associated with S x , S y , whereas Hare associated with s x ,s y ). 

It will be noted that the values of s x and S x for plate 6 are exceptionally large 
(9.1 and 7.8/im, respectively). Referring to the plotted results in Figure 12, we see the 
presence of two pronounced jumps in the trends as indicated by the heavy arrows, lo make 
sure that these jumps were actually real and not the result of a measuring or recording 
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Table 1. RMS values of x and y drift's before and after regression. 


Plate 

RMS Values Before 
Regression (V-m) 

s x s y 

RMS Values After 
Regression (pim) 

S, s y 

1 

3.4 

2.5 

2.4 

1.8 

2 

3.5 

2.3 

3.4 

1.8 

3 

4.0 

1.7 

3.7 

1.9 

4 

5.3 

2.3 

4.3 

1.1 

5 

3.2 

2.1 

3.3 

1.9 

6 

9.1 . 

2.3 

7.8 

2.2 


£6x^1 

i 

s 

S. = 

[ClSxj -Sx-,) 3 ! 

11 

X 

8 

poyfl 

1 

2 

s = 

fS(6y i -6 y ;f'l 

l n 


8 1 


6 x-,/ 6y - observed values 

J J 

6 y j = values computed from regression 
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ifferent comparator by different operator. 




blunder, we had plate 6 remeasured by another operator on another comparator (namely, 
the DBA multiiaterative comparator). Comparative results presented in Figure 13 show 
good agreement between the two sets of measurements and verify that the jumps are 
indeed real. Reviewing the plots for the other plates, we find indications of similar 
jumps occurring on Plates 2, 3, 4 and 5 at the points indicated by the arrows in Figures 
8,9, 10 and 1 1 . These suggest the presence of stiction in the drive train. Stiction being 
the occasional slight binding of the gears followed by sudden release. Whether or not this 
is the correct explanation for the abrupt changes remains to be determined. For convenience, 
we shall nonetheless refer to the abrupt changes as stiction jumps . 

It seemed that the generally mediocre fits obtained from the regressions performed on 
6x could well be attributable to the compromising effect of stiction jumps. To test this 
hypothesis we revised the regression model so that the single zero set term a Q was replaced 

by the expression: 

Si a 0 + Ca a ° + £3 a ° 

in which, 

£1 = 1 for Tj < Tj_ 

= 0 for r t > T 1 

£ 8 = 1 for T x < < T 2 

= 0 otherwise 

= 1 f o r T 4 >T 2 

= 0 otherwise. 

This artifice makes it possible to reinitialize the zero set term up to two times per plate. 

When no reinitializations are desired, we need merely set T x and T 2 to values greater than 
the maximum value of n (say to 180 seconds). The values of T x and T 2 adopted in the 
revised regressions for 6x are as indicated in the following table. 
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Plate 

Number of 
Stiction Jumps 

Ta 

t 2 

Degrees of Freedom 
For Regression 

1 

0 

180, 

180 

8 

2 

1 

7.5 

180 

7 

3 

1 

22.5 

180 

7 

4 

1 

37.5 

180 

7 

5 

1 

127.5 

180 

7 • 

6 

2 

67.5 

142.5 

6 


The results of the revised regressions are plotted in Figures 14 through 19. The 
results for Plate 1 are unchanged, which is as It should be, since no stiction jumps were 
exercised. Results for the remaining plates represent substantial improvements over those 
from the original regressions. The fitted functions now very closely follow the observed 
values. This is verified by Table 2 in which Table 1 is extended to reflect the results of 
the revised regressions. The grand rms value for the fit in 6x for all six plates is 2.1 ftm, 
a value only insignificantly greater than the grand rms of 1 .8/im for y. 

The amplitude of the periodic error from the revised regressions averages about 0.6 
second of arc in right ascension and about 0.4 second of arc in declination . This is 
appreciably lower than the values considered typical (i.e. 2 to 3 arc seconds in right 
ascension) in Reference 2. However, specific results from the MOTS at Fort Myers wpre 
not reported in Reference 2. 

1.6 GENERAL CONSIDERATIONS 

When stiction jumps are duly taken into account whenever they occur, excellent 
fits are obtained from the simple regression model that was adopted. This demonstrates 
the effectiveness of the rotating diffraction grating as a means for monitoring the stability 
of the MOTS camera throughout routine operations. 

Because of the exploratory nature of our investigation, we took the liberty of 
ignoring a few fine points that should be implemented when the diffraction data are 
actually to be used to generate corrections applicable to satellite observations. Because 
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Table 2. Extension of Tabie 1 to include results of revised regressions. 



6x , 6/, = observed values 

3 3 

6x', 6y' = values computed from original regression 

3 3 

5 X " f 5y" = values computed from revised regression 

3 3 

f = degrees of freedom 

= 8, 7 or 6 depending on number of stiction jumps 
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the length of each exposure through the grating may be a significant fraction of the 7 90 
second period of the drive, the integrated effect of the drift over exposure increments 
should be recognized as being what is actually observed. Accordingly, if 2A denotes 
the length of each exposure through the grating, the expression 


. 2rr 277 

a 2 sin 90 T - ! +a 3 cos 90 r i 


in the formula for 6x t , for example, should more properly be replaced by the expression 


, f T ‘ +A r . 2JT 2lT . 

zk Jr.-A ta2!m 90 t+ ° 3COS 90 ,]d ' 


However, this reduces to the form 


. . 277 , , 277 

a s sm 90 r i + ° 3 cos 90 T i 


in which 


°2 ^ 


K a. 


where 


K = (sin % a)/(|( A) 


Thus,while a^ and a^ are the values actually obtained from the least squares regression, 
the values to be used in the correction formulas should be a 2 ,a 3 which will, of course, 
be slightly larger. For 2A = 15 seconds (the value used in our experimental investigation), 
the value of K is (sin T7 /6)/ (tt/ 6) = 0.956, which means that the amplitudes obtained from 
the regressions should be increased by a factor of 1/.956 = 1 .044. 
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The second fine point to be considered in operational applications of the 
diffraction grating concerns effects induced by changes in atmospheric refraction 
resulting from the gradual changes in zenith distance occurring throughout the operation. 
Even if the MOTS were perfectly stable in the right ascension-declination frame, the 
gradual change in refraction with change in zenith distance would be' manifested as a 
secular drift if appropriate corrections were not applied. In Reference 3 we showed 
that the correction for this effect is given by: 

d(AC) =4.5x10 t t sin cos sec 2 Cj 
(arc sec) 

in which, for the present application, 

angular correction to be applied zenith distance of mean of pair of 
diffraction images 

time of exposure of i th set of diffraction images relative to mean 
time of total exposure 

azimuth of star at t t 

zenith distance of star at t t 

latitude of station . 

By means of standard astronomical formulas, this correction for zenith distance can be 
propagated into right ascension and declination and thence into x,y plate coordinates. 

It is the corrected x,y coordinates that properly should be employed in the diffraction 
grating regressions. 

1.7 CONCLUSIONS AND RECOMMENDATIONS 

The diffraction grating apparatus evolved from this study has proven to provide a 
simple, inexpensive, effective,and operationally unobtrusive means for monitoring the 
stability of the MiOTS camera to a precision of a few tenths of a second of arc. 


d(AC) = 



$ 
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We recommend that NASA employ the breadboard apparatus developed by DBA 
in further tests on other MOTS cameras, particularly those that are known to be affected 
by I arger drift errors than the MOTS at Ft. Myers. Should these tests confirm the promise 
of the approach, NASA should undertake the procurement of operationally optimized 
units to be used routinely in future operations at all MOTS stations. Such units should 
preferably be controlled by an automatic programmer with a suitable range of selectable 
exposure rates. 

While the diffraction grating method provides a way of correcting for drift in future 
MOTS operations, the question naturally arises as to whether anything can be done about 
the many hundreds of plates gathered on past operations, particularly on GEOS I and 
GEOS II. We believe that quite possibly something effective can be done. Because a 
formal standardized procedure was followed in the exposure of GEOS plates, a fairly 
high level of repeatibility of drive error may well exist. In any event, whether or not 
this is the case can be determined experimentally by employing the grating in repeated 
trials simulating GEOS operations. If acceptably repetitive results are obtained for a 
given camera, a pooled result can be used to derive corrections to be applied to directions 
obtained from previous plates taken by the camera . The validity of such corrections can be 

tested by determining whether or not improved residuals are generally obtained from short 
arc reductions based on revised directions. 
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SECTION 2 


MINITRACK DATA REDUCTION 
2.1 INTRODUCTION 

This section covers work done by DBA Systems, Inc., over a 2-year period. 

Work for the contract was initially being performed by the DBA Florida Office. 

In order to provide better liaison betweem NASA and DBA, it was decided to transfer 
the project to the DBA Washington Office. This was done nine months after the start 
of the contract. The section dealing with the Minitrack geometrical error model, 
contributed by James B. Willmann, was written prior to the transfer of the project to 
Washington. Several thousand NAP (Network Analysis Program- II) control cards for 
processing optical orbits had also been prepared, but this work became largely redundant 
after the development of a new optical preprocessing program. 

After the project was transferred to Washington, work was concentrated on cali- 
brating Minitrack using short optically determined reference arcs. This was done in 
order to demonstrate the feasibility of using the NAP program for the purpose of cali- 
brating Minitrack. At this time too, the magnitude of the Doppler effect was being 
investigated, and the relevant equations were implemented into the NAP program. 

Concurrently, the new optical preprocessing program, already mentioned, was 
developed. The reason for writing a new optical preprocessing program was that the 
old program could not process data for SAO stations and required an enormous amount 
of data preparation. The new program in addition would punch out most of the NAP 
control cards. 

Twelve months after the start of the contract, work began on the calibration of 
Minitrack using Minitrack data only. Initially, little progress was made. At this time 
the Minitrack extract and sort programs were developed. The Minitrack preprocessor 
was modified to punch out most of the NAP control cards. After three months, it was 
found that the NAP program had been requested to compute and apply the effect 
of the earth's precession and nutation. Unfortunately, an error was known 
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to exist in these computations. After the NAP control cards were changed, so as not 
to request these computations, progress was again made. (To the order of the accuracy 
of the Minitrack system, it would not appear to be necessary to consider precession 
and nutation. However, because of the surprisingly good agreement between orbits 
determined from Minitrack and optical data, this point of view may have to be 
modified). 

Definite but slow progress was now being made. A post NAP program — later 
extensively modified — was written to summarize the results of the NAP program. This 
also involved modifications to NAP subroutines RESID and FINALP. 

About 17 months after the start of the contract, the contract Technical Monitor 
(W. M. Rice) proposed an analysis of the NAP program by Boole & Babbage, Inc., 
to determine the feasibility (within the time scale of the contract) of increasing the 
computational speed of NAP and hence the overall rate of progress. The analysis 
showed that the majority of the computing time was being spent in a few relatively 
short subroutines. A new algorithm was then developed for the computation of spherical 
harmonics and this was implemented into NAP. The modifications to the NAP program 
doubled its overall computational speed. 

Advantage was now taken of the faster NAP program, the length of arc con- 
sidered was increased from 2.5 to 5 days; however, this led to further problems. 

A meeting was called by E. P. Damon of the Computer Systems Branch to call attention 
to the amount of computer disk space being required by the NAP program. (As the 
length of arc had been doubled, so had the amount of data being stored on disk 
memory). We had, however, become aware of the problem before the above 
mentioned meeting because of the frequently aborted runs due to non-availability 
of enough disk space. To overcome this problem, the special subroutines PRTIAL and 
RESID were written and the NAP disk space requirements were substantially reduced 
The post-NAP program was modified to be consistent with the latest NAP program . 
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With the increasing amount of data being processed on each run, the number 
of NAP control cards became very large and it was decided to write them on tape. 

To handle modifications to individual NAP "cards" the pre-NAP card updater 
program was written. 

In checking out the new special subroutine PRTIAL a discrepancy was discerned 
in the printed out time. This was traced back to the Minitrack preprocessing program, 
which was corrected. All Minitrack data used in the data processing now had to be 
reprocessed . 

At about 21 months from the start of the contract it was decided to investigate 
the puzzlingly slow rate of convergence of the solutions computed by NAP. The 
investigation revealed that the NAP matrix inversion was not being handled as 
described in the NAP documentation but in an iterative manner, which was extremely 
time consuming for the long arcs being considered. The complete equations as given 
in the existing NAP documentation were incorporated into NAP by R. Garza-Robles 
of Goddard Space Flight Center (Code 551). The rate of convergence was at least 
doubled . 

With the help of the support programs already described and with the modified 
NAP program a four-arc run was now made on the 360/95. The total time span covered 
was 25 days and the computer time taken was 20 minutes. (At the start of this program, 
one 2.5 day arc was processed in 24 minutes and then the convergence rate of the 
solution was half of its present rate). "Optical" reference arcs were computed for 
the same time spans and a comparison of results was made. Advantage was taken of 
the current DBA contract (NAS5-11730) to develop an ionospheric error model for 
incorporation in the NAP program and one run was made using ionospheric corrections. 

It should also be mentioned that the NAP-II program was successful ly overlaid 
reducing its core requirements from 700,000 to 500,000 bytes. 
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This section is devoted to discussion of the reduction of Minitrack data using 
the NAP-II program and auxiliary programs written specifically for this task. Included 
in this section are discussions of the Minitrack error model and the Doppler effect 
and wave-propagation time delay as applied to predicted Minitrack measurements. 

The second and third parts of this section discuss the handling and preprocessing 
of the Minitrack and Optical data respectively. The fourth part is devoted to the 
reduction of the data using the NAP-II program. 

Results for a reduction effort are given in the fifth part. Results are given for 
the station at Ft. Myers and then a multi-station network which includes all of the 
Minitrack stations. 

The last part of this section is a recommended procedure for reducing Minitrack 
data on a production basis. 

As of special note: Appendix A-4 should be continuously referenced in the 
data processing cycle, as it contains restrictions on the programs and assumptions 
pertaining to the data. 
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2.1.1 The Mini track Geometrical Error Model 

This discussion will consider some of the general geometric properties of the 
Minitrack system. As with any interferometer system, the basic measurement mode 
is that of the phase path difference between the target and two antennas located 
some distance apart. The line between the antennas is referred to as the base line. 
The geometry for two cases is depicted below; 'A' for a target at a close range, 
and 'B' for a target at infinity. 



A B 

Figure 1 


In Figure (1), A and B, the antennas are at the end of the base line a. The 
measurement made by the system can be interpreted as being proportional to the 
distance P. In both cases, P is defined as 

P — Rg — Rl . 

Generally, we say that the system measures the cosine of the angle 9 

since 

P 

cos 9 ~ — . 
a 
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For the case of a point at infinity, the above expression is exact. For the 
case depicted in Figure 1A, it is only approximately true and cos 0 depends upon 
R as well as P and a. 

For t!.-3 Minitrack system, there are fcur antennas located at the ends cf 
two perpendicular base lines. These base lines are usually located to intersect 
each other at their centers. However, since exact location of the phase centers 
of the antennas is not possible, there will be errors in the resulting data if they are 
interpreted as a true measurement of the cosine of the direction to the target. 
Therefore, the discussion which follows will investigate two areas. First, we will 
derive the relationship between the angle 0 and the measurement P. Then, we 
will investigate errors in the antenna locations and how they effect the interpretation 
of the data. 

To derive the desired expression, we can refer to Figure 1A. From this 
figure, the following relationships can be written. 

P = R 2 - R, 2 

Ri 2 = R 2 + (p ■) - R a cos 0 

R 2 2 = R 2 + j + R a cos 0 . 

We will redefine the measurement as 



Then t 0 approximates the cosine of 0. Proceeding with the algebra, give 
P 2 + 2PR X + R x 2 = R 2 2 
P 3 + 2PRi - 2Ra cos 0 = 0 

(P 2 - 2Ra cos 0) 2 = 4P 2 R x 2 = 4P 3 (R 2 + ^ - R a cos 0) 

F' 1 + 4R 2 a 2 cos 2 0 = 4P 2 R 2 + P 2 a r 

solving for cos 0 give 

— — — — — 

cos 0 = lo y 1 1 + ^ 2 - (1 - -to 2 ) 
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the desired result. Note that as R-*«> then -to-*cos 0. Expanding gives 

2 

COS 0 — -t-o 0^2" (1 “ ) • • • 

The maximum error in assuming cos 0 = t 0 occurs when t Q is approximately 
l/v/3 - . Th is leads to a maximum error of about 

a . ~ a 3 

cos 0 - -Cq - 20R 3 * 

For this analysis, we will consider the three dimensional case pictured in 
Figure (2). The antennas are located along the x and y base lines which are 
approximately a and b in length. We will then relate the path difference measure- 
ments to the desired cosines as a function of a, b and the errors in antenna location. 
This is pictured below. 



Figure 2 
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We have also included a figure which defines fhe measurements of azimuth 
and elevation. First a few definitions. 

I = cos a 
= cos jS 

n = cos y - \! 1 - cos 3 a - cos 2 /3 = / 1 -l 3 - m 2 

x = R • -L = R • cos A cos E 

y = R • m = R • sin A cos E 

z = R • n = R • sin E . 

Thus, 

tan A = m /i ; cos E = / l 2 + m 2 . 


The measurements made by the system are given by 

P = R 3 - Ri = V (x + ^ - x 3 ) 3 + (y - y 3 ) 3 + (z, - z 3 ) 2 

- /(x - 1 - x x ) 3 + (y - y x ) 2 + (z - z x ) 2 

Q = R 4 - R x = n/(x - x 4 ) 2 + (y + j - y A f + (z - z 4 ) 3 

- '/(x - x 2 ) 3 + (y - | - y 2 ) 3 + (z - z 3 ) 3 . 


In the above, we have chosen the coordinates of the antennas to be 

1 x x + a/2, y x , z x 

2 x 2 , y 2 + b/2, z 2 

3 x 3 - a/2, y 3 , z 3 

4 x 4 , y 4 - b/2, z 4 . 

In the case where there are no errors, then x x , y x , ..., y 4 , z 4 would all 

be zero. In that case, if P 0 is the value of P when all survey errors are zero, then 

P 0 = \/(x + a/2) 3 + y 3 + z 3 - %/(x - a/2) 3 + y 2 + z 3 

= A 3 + ax + a s /4 - /R 3 - ax + a 2 /4 
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and 


Qo = n/R 2 + by + b 2 /4 - \/R 2 - by + b 3 /4 . 

The solution to these two^equations is the same as that presented in Section 
II. That is 


, _ Po 
1 ~ \ 


c4) ■ 

Po , a s 

- m 

£(> 

P 0 2 \ 

~’3 r ) 

Qo 

m = nr) 

/ 1 4R a 

>• -%f> ■ 

Qo , b 2 

"T + 8F 

Qo 1 
"T \ 

' Qo 2 

} ~P" 


However, of more importance is the effect on l and m when the errors in 
antenna location are not zero. To explore this, we can expand the equations for 
P and Q about x x ... z 4 =0. This gives 


but 


thus. 


P = Po - 


(x + a/2) X 3 + yy 3 + zz 3 (x - a/2) X x + yy x + -z.-z.-y 


R3 Ri 

For purposes of this expansion, we can let R 3 = R x = R, giving 

p = p ° + ^ ( x i " x a) + 1 (yi - Vs) + f ( z i " zs) ’ w ( x i + x s) 


a x y 

1 " R ' m " R ' 


R 


R 


2R 


z 

n = R 


p = Po + - x 3 ) + m(y x - y 3 ) + - z s) - ^ ( x i + x s) • 


The actual measurement is given by 


= a ' 


and from before 


~ £(1 - -t 3 ) « £ . 

Therefore, we finally have 

. - » . K - x 3) _ (yi - ys) _ ( z i - z 3 ) L ( x i + x 3 ) 

- -Co - m 0 - n 0 + *5 

a a a zK 
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Proceeding in the same manner, we have for m 


_ _ _ , (x 2 " X J 

m — mo — -Co — — r 


.. (y 2 - yJ _ (z 3 - Z 4 ) J y 2 + y 4 

~ m ° — E n ° — E 2R — 


Normally, the last term of both expressions can be neglected. Alternately, 
if a and b are chosen so that x x + x 3 and y 2 + y 4 are zero, then it can be dropped 
from consideration. 


For close-in ranges when the approximation 
a 

cannot be made, then the value of t and m above can be used to compute the exact 
cosine value through the solution of 


1 = + 


a 2 P 2 


o n . Po 2 \ 

a ( ‘ IT) . 
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2.1.2 The Doppler Effect and Wave-Propagation Time-Delay as Applied to 
Predicted Minitrack Measurements 

Consider a satellite at position r(t) at time t transmitting an unmodulated 
signal at a constant frequency v . The signal is received by two MINITRACK antennas 
situated at ^i and - ^ ? , respectively, where a is the baseline length and i_ is a 
unit vector in the x-direction. The phase difference between two signals arriving at 
the two antennas at the same time t is measured. If the two signals left the satellite at 
times t A and t B , respectively, then the phase-difference between the two signals is 
2 77 1/ (t A — t B ) . The actual measurement, A^) , is the phase-difference divided by 2 77 . 
Hence, 

At p = v (t A - t s ) . (1) 


or equivalently. 


A <p 


_ c (t A - t 0 ) 


( 2 ) 


where c is the velocity of light and X the wavelength . We must also have. 


llA - |il = C ( f " f A) ' 


(3) 

and J_r@ + ^ j_| = c(t - t B ) , 


(4) 

where, in general, r t = r (tj 



Writing, 



r = - *0 »• 


(5) 

we have by the mean-value theorem 



U = U ~ r r A + £ t 2 r*„ 

/ 

(6) 


where t M is a time between t A and t s . 
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Writing 


( 7 ) 


r A =|la| Qnd Pa = J>Aa , 

we obtain from (6), 

ij + lai 


Since p A 2 = 

1 / Pa *Ja = ? A / and _i • p A = I A , 


the direction 

cosine at time t A , we obtain from (8), 


1 

+ “ r *{ ,+ ( r )[ l ‘- 2 'o c + '» 3 



+ (r 4 ) [*i- '« c + *'.•?-] 

| i f m 

where 


(10) 



Next we make some assumptions regarding the magnitude of the terms in ( 9 ). 

We assume that, 

(la I £ 1 (II) 

(Th is will be justified later on when we shall show that l 0 is a first 
approximation to the direction cosine l A ) 

Write, 7 = a 3 , ( 12 ) 

'A 

and assume that a s < 10“ 3 , i.e. that the satellite is always more than 1000 
baseline lengths from the station. 

Write = u 4 ( 13 ) 

c 

and assume that | | /c: < 10 -4 , i.e. that the satellite velocity is less than 
10 -4 times the velocity of light. Since |? A | =Jr A | it follows that J U4 J < IQ " 4 . 
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( 14 ) 


Next assume that, |r| < 1CT 7 
Hence by the assumption made In (12) , 

P liJ < w 10 O 5 ) 

With the above assumptions equation (9) may be rewritten as, 

|r, + tai_| = r,{l +a s [l,-2l„ u 4 + a 3 + 0 (10-^)| } * (16) 

4 

where 0 (10~ n ) denotes terms of order 10 -n or less. 

Equation (16) may be rewritten in the form: 

|r s + 5 aj_ | = r,{l +° 3 l.+ jc.f I, 3 + a 3 [l a 3 (1 - l,=) - 2 1 0 u 4 + 0 (1CT 7 ) 

= r,(l+|a 3 l,){l+a 3 [l a 3 (1 - l A 2 ) - 2 l 0 u 4 +0(10 -7 ) 1 -c 3 l,+0(10- 6 ) 

= r l (l+ 50 3 l J ){l+^ o 3 ^a 3 (l-l, a )-2l 0 u 4 +0(10-’) l-a 3 l A } 

= r.{l+5a 3 l,+ja 3 ^a 3 (l-l,=)(l-2a 3 l,)-2l o u 1 +0(10-') ) (16)' 

We obtain a similar expression for [r x - i a j |with a 3 replaced by -a 3 and u 4 = 0. 
Hence: 

|r A " 5 ° 1 l = r »{’ - ^ °3 + 5 a s “sf 1 - | « a )(' + 5°3 l » 

+ 0(10-’)]} (17) 

Hence from (3), (4), (5), (16)) and (17), 


c t - r A a 3 |l A - 1 a 3 3 (l-l A 3 ) l A - l 0 u 4 + 0 (lO -7 )} (18) 
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It is easy to show that ^ a 3 2 (1 - l A 2 ) l A = 0(10~ 7 ) . Hence from (10), (12), 

o 

and (18) 

lo = Ia - >0^4 + 0(10- 7 ). 

Whence, by (13) 

I, = l 0 (l +-f) + 0(10- 7 ) (19) 

Hence we have shown that l 0 is a first approximation to the direction cosine l A . 
From (2), (5), and (10) \ 

| = Ms? , (20) 

0 a 

and hence 

'* ■ f - [*<'♦£> 

From expression (21) it is easy to see that the term in (r/c) represents the 
Doppler effect. We shall show that expression (21) may be further simplified. 

From (3) and (17) we have 

c (t-t A ) = r A 1 +0(10- 3 ) (22) 

Hence denoting the direction cosine at time t by I, we have by the mean 
value theorem 

where l N is the direction cosine at some time t N between t and t A . Writing 
p = r/r we obtain 

• • • 

r p - _r - r p 

. = p x (r x p ) (24) 

since p 3 = 1 and p • jr = ? . Differentiating (24) gives, 

r p = p x (r x p ) + px (r'x p ) + px (rx p ) - r p (25) 
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Since I = p • i If follows that. 


1 s |E|- 

1 =|i 

) | and 1 = j p 

1 

By (24), 


r |p| 5 |f_| 


From (25) and (27) 



r 2 

p | = 

r | 2 + r 1 ¥ j + 

i\ a + f |i 


< 

3 | r 1 2 + r | 

1 — 1 1 

r [ 


From (26) and (27) we have 



c c 


= 0 ( 10 “ 4 ) 


by our previous assumptions. From (26) and (28) we obtain 


1 

2 


SE £ 


o i r I 3 + r | r 1 

J c 2 c 2 


= 0 ( 10 “ 7 ) , 


by the assumptions made in (13) and (14) above. 


We then have from (23), (29), and (30), 

I = l A +— |‘ A + 0(1 0 -7 ) + ^4 0 (10- 7 ) 

c r N 

By the mean value theorem 

r N - r A + On *a) r x 

But t N -t A ^ t - t A =^ A 1 +0 (10“ 3 )] , Hence 

| 1*n - t A ||f K | ^ r A |r K |/<= l 1 +0 ( 10_3 )J 

= r A 0 (10 -4 ), where 

r N > r A - r A 0(10 -4 ), and (r A /r N ) = 1+0 (10' 4 ) 
Hence I = l A + ^- A l A +0 (10 7 ) 
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( 33 ) 


But r A h X A “ r A I A / 

where \ - _r a * _[ 

Hence I = l A + - l A (A + 0 (10" 7 ) (34) 

From (19) we obtain 

, o = | A (1-1 A ) + 0 (10— ) , 

Hence from the above and (36) 

1 = 1 + ^ + 0 (icr 7 ) (35) 

0 c 


By the mean value theorem and (22) 


x 

c 



i + o (icr 3 ) 


= - A + 0 (icr 7 ), by (14) 
c 

Hence, £ 

i 8 = I - c- +o (icr 7 ) . 

Finally from the above and (20), 


where x is the x-component of the 
velocity at time t. 


= | - - +0 ( 10 ~ 7 ) 

a c 


(38) 


The following assumptions were made in deriving (38) 


(a/r) < Kr 3 , ([r|/c) < KT* , 

and r j r f /c 3 < 10 -7 . If the only forces acting on the satellite are due to 
the Earth's gravity field then |r_| = gR E 3 /R 3 , where g is the acceleration due 
to gravity at the Earth's surface, R E is the radius of the Earth and R is the radius 
vector from the center of the Earth to the satel'ite. Since r ^ 2R it then follows 
that r jrj/c 3 S [2g Re/c 2 |r e /r| which is always less than 10~ 7 . 
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For the m direction cosine a formula similar to (38) holds true: 


X A'4 


m 


= m _JL + 0(10 -7 ) , 

c 


( 39 ) 


where in this case 

m 


represents the phase difference of the North-South antenna pair. 
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2.2 MINITRACK DATA 


2.2.1 Standard Minitrack Data Format 

The following is the Standard Minitrack data message that is stored on magnetic 
tapes by satellire and by station. The message appears on teletype output as in the 
following sample . 


Satellite ID 
and date 


Data 


Satellite ID 
Terminating 
Message 


& 64 064 01 1 6 
4350.26391 14 
1456.307 175 0 
1662.3591 152 
1861.40? 1559 
2061.4 74 1972 
2263.5342388 
24 63.591 1798 
2664.65 022 08 
2863.7 162 633 
3064.7882057 
32 65 .55 724 83 
3467.94829 18 
3668.0162345 
3867. 102277 1 
4 068. 18922 00 
4269.233 1634 
4469.372 1036 
4 670.4 74 15 14 
4870.5701966 
507 1.669 1395 
527 1.7781843 
54 72.889 1239 
5672.934 175 1 
5874 . 1 14 02 18 
0074 .234 155 3 
0274 .3590109 
0475 .47505 68 
067 6. 62 00003 
0879.7680522 
1073.3350954 
1230. 102 0389 

64 0 64 01 
03/ 132 1Z JAN 


90103 

.3580.2639 1 14 
.4503.323 1736 
.4505 .3581 199 
.4 5 04 . 4 2 4 15 9 7 
.4505 .4672022 
.4506.539 1436 
.4507.59 12840 
.4507.649 1255 
.4505.72 12679 
.4508.7992 104 
.4507.8672528 
.4507.9431960 
.4508.0281380 
.45 09. 1 1623 12 
.4507.201 1247 
.45 10.2992684 
.45 1 1 .387 1 1 15 
.45 1 1.432 15 60 
.45 1 1.530 1 993 
.45 13. 634 145 0 
.45 12.734 1909 
.45 14.896135 1 
.45 15 . 003080 1 
.45 16. 1 190238 
,4 615.244 0702 
.4 616.3 74 0175 
.4619.5000613 
.4 619.6460072 
.4 620. 7230527 
.4 620.82300 14 
.4 62 1.098042 6 


L W N K 


. 1230.2639 1 14 . 
. 12 17.3 121830. 
. 1222.5761243. 

. 1217.422 1636. 

. 1220.4761060. 
. 1222.5501463. 
. 1227 .6082386. 
. 1228.6522305. 
. 1232.75 127 17. 
. 125 6.3062 143 . 
. 1233.5842574. 
. 1235.9612005. 
. 124 1.037243 1. 
. 1247. 1242870. 
. 1249.2 102297. 
. 125 1.2932731. 
.1256.5921171. 
. 1259.492 1607. 
. 1261.5871055. 
. 12 65 . 693 15 03 . 
. 12 65 . 795 1947. 
. 127 1.9 15 1597. 
. 1275.023 1826. 
. 1278. 1570252. 
. 128 1 .27007-19 . 
. 1233.399021 1. 
. 1287.5 170665 . 
. 1290. 633012 1. 
. 1294 .79505 57 . 
. 1293.93 10035 . 
. 1200.0230478. 


00330.2659 1 14 
00554 .325 1831 
00355.377 1280 
00339 .4252685 
00342 .48 12033 
0C347.54 1 15 19 
00343.6062920 
00350.672 1342 
0035 1 . 7372762 
00354 .3202 186 
003 5 8. 8842 604 
00359.959203 1 
00366. 04 S24 59 
003 65 . 13 129 04 
00366.2 152342 
0037 1 .3 13 1768 
00370.403 1207 
00376.494 1661 
00376. 59310S7 
00379.694 1544 
00353.3001990 
00335.9 15 1439 
00337.0161901 
00390. 169 1333 
00393 .2360777 
00397.3930246 
00300.4330695 
U03G2 . 654 0176 
00305 .7900650 
005 08 .9 62 0 12 6 
005 09.05 005 07 


.2 15.2639 1 14 . 
.2 15 .3391904 . 
.2 15 .375 1324 . 
.215 .4321729. 
.2 15 .4831 142 . 
.2 15 .5472546. 
.2 15 .6102969. 
.2 15 .6792535 . 
.2 15.7452302. 
.215.8192239. 
.2 15.3961652. 
.2 15 .9802073. 
.215.05825 17. 
.2 15. 14 12946. 
.2 15 .2252584 . 
.2 15 .3 05 1830. 
.2 15.4 12 1255. 
.2 15.5021704 . 
.215.603 1146. 
.215.7051598. 
.215.8281045. 
.215.9510481. 
.215.0611942. 
.2 15 . 1500594 . 
.2 15 .285 0c35 . 
.215 .407 0299 . 
.2 15 .5330745 . 
.2 15 .647023 1 . 
.2 15 .3420678. 
.2 15.9890143. 
.215.0990595. 


The format for the message is as follows: 
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MESSAGE FORMAT 


Characters) 


Contents 


1 , 2 
3, 4 
5 

6 to 8 

9 

10 to 12 
13 
14, 15 
16, 17 
18 

19 to 21 
22 

23 to 25 
26 
27, 28 
29, 30 
31 

32 to 34 
35 

36 to 38 
39 

40 to 42 
43, 44 
45 

46 to 48 


second (time) of frame start 

hundreds and tens digits of east-west medium phase 

period "."separator 

first east-west fine phase 

signal strength indicator (AGC) 

first north-south fine phase 

period " . " separator 

minute of frame start 

hundreds and tens digits of east-west coarse phase 

period separator 

second east-west fine phase 

signal strength indicator (AGC) 

second north-south fine phase 

period "." separator 

hour of frame start 

hundreds and tens digits of north-south medium phase 

period " separator 

third east-west fine phase 

signal strength indicator (AGC) 

third north-south fine phase 

period " . " separator 

day of year for frame start 

hundreds and tens digits of north-south coarse phase 
period " . " separator 
fourth east-west fine phase 
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MESSAGE FORMAT (Cont'd) 


Character(s) 


49 

50 to 52 

53 

54 


55 , 56 


57 

58 to 60 

61 

62 to 64 
65 


Contents 

signal strength indicator (AGC) 
fourth north-sound fine phase 
period " separator 
equatorial/polar antenna indicator 

1 = equatorial 

2 = polar 
station number 
03 = FTMYRS 

05 = QUITOE 

06 = LIMAPU 
08 = SNTAGO 
12 = NEWFLD 
period " . " separator 
fifth east-west fine phase 
signal strength indicator (AGC) 
fifth north-south fine phase 1 
period "."separator 


15= WNKFLD 
16= JOBURG 
19= ALASKA 
21= ORORAL 
23= MADGAR 


01 = BPOINT 

13 = COLEGE 

17 = MOJAVE 

14 = GRDFKS 

18 = WOOMER 


2.2.2 Extracting and Sorting Minitrack Data 

The GEOS-A Minitrack data extracting program and the s.orting program are 
special purpose programs each with a single function. The extracting program reads 
a Master Minitrack data storage tape and extracts only data taken on the GEOS-A 
satellite. This data is copied to another magnetic tape in the same format (same as 
given in the example in 2.2.1). This results in a data tape of GEOS-A Minitrack 
data grouped by station in chronological order. The FORTRAN listing of the 
extracting program is given in Appendix A-2.1 . The requirements for the program 
operation are: 


53 



Function 

Unit 

Description 

Card Input 

Card Reader 

Program Source Deck 

Tape Input 

Data Set 9 
(FT09F001) 

Mini track messages for several satellites 

Tape Output 

Data Set 1 1 
(FT11F001) 

Minitrack messages for GEOS-A 
(Satellite ID 65891) 

Printed Output 

Data Set 6 

Total number of messages 
Number of GEOS-A messages 
First 100 messages from output tape 

The tape output (Data Set 11) from 

the extract program is the input data tape 


for the sort program. This program sorts messages in time-sequence. It also eliminates 
duplicate messages. Messages for the same station and antenna configuration are 
assumed to be duplicates if they commence within 30 seconds of each other . If two 
duplicate messages are of unequal length, then the shorter message is discarded. If 
they are of the same length, then the last received message is discarded. 

The length of a message is determined by the number of good records it contains. 
The characters in each record are checked for numeric characters and periods in the 
appropriate places. A record with an error is discarded. The first "proper" record 
of a message is the calibration line. If this contains an error, the whole message is 
discarded. It has been found that in some messages the calibration line and sub- 
sequent records are repeated within the same message. For this reason, every record 
is compared with the calibration line. If a record is found to be identical to the 
calibration line, then this is regarded as the first "proper" record of the message. 

All previous "proper" records are ignored. 

Messages containing fewer than 5 records are discarded. The output tape from 
the sort program is in the same format as the input tape. 
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The FORTRAN listing of the sort program is given in Appendix A-2.2.- 


The requirements for running the sort program are as follows. 
Function Unit Description 


Card Input 
Data Cards 


Card Reader Source Program 

Data Set 5 17 Cards Required (See Appendix A-2.2) 

Cards 1-16 (One card for each Station) 


Card Column 

Word 

Definition 

Format 

1-5 

KSTA 

Station ID 

15 

6-10 

11-16 

STATIO 

Blank 

Station Name 

A6 

Card 17 

Card Column 

Word 

Definition 

Format 

1-8 

INTAPE 

Input tape No 

. A8 

9-13 

INFILE 

File No. 

15 

14-21 

NDTAPE 

Output tape No.A8 

22-26 

NDFILE 

File No . 

.15 


Tape Input 


Tape Output 
Printer Output 


Data Set 9 
(FT09F001) 

Data Set 10 
(FT10F001) 

Data Set 6 


Input tape containing Minitrack messages 
(Output tape from Extract program) 

Output tape containing sorted Minitrack 
messages 

A list of message numbers (as sequenced on 
input tape) of rejected messages together 
with reasons for rejection 

A list of output information giving the 
following. (The actual data message is not 
printed out): 

Sequence number on output tape 
Sequence number on input tape 
Station ID (KSTA) 

Station Name (STATIO) 

Day Number (1 to 366) 

Month 

Day of Month (not correct for leap: year) 
Seconds of day ~ - 

Hours | 

Minutes >- (of day) 

Seconds I 

Number of lines in message . w ,. ■ 
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2.2.3 7 he Mini track Preprocessor 

The purpose of the Mini track preprocessor is to make known corrections to the 
Minitrack data and output a magnetic tape containing corrected data in a format 
acceptable to the NAP-II program. Basically, the preprocessor reads the raw phase 
data message, computes phase differences, makes corrections and then converts the 
corrected phase differences to direction cosines. This process is explained in detail 
in references 1 and 2. 

The Mini track preprocessor used was based on that developed by the Goddard 
Space Flight Center Network Computation Section (references 1 and 2). The program 
was extensively modified by W. M. Rice (GSFC, Code 551). A few minor modifi- 
cations was then made by DBA Systems, Inc. under this contract. A complete listing 
of the preprocessor program is given in Appendix A-2.3. Sample JCL for running 
the program is also given. 

Changes Made by W. M. Rice . A Minitrack message contains, in 
general, 155 measurements of phase differences between two antenna 
pairs. The preprocessor converts these to 155 pairs of direction cosines. 

A quadratic fit is then made to each set of 155 direction cosines. The 
resultant output is a single pair of direction cosines (the fitted mid- 
point values). The program is made more complicated by the fact 
that the ambiguity of the direction cosines must be resolved using 
additional information in the Minitrack message (Reference 2). 

Flowever, this is only done to the fitted midpoint values . In the 
Rice version of the program the ambiguity is resolved for all 
direction cosines. The output consists of all direction cosines and 
is in a format that can be used as an input to the NAP program. 

The time associated with each measurement is given at the instant 
that the signal being analyzed arrived at the relevant antenna pair. 
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Changes Ma de by DBA Systems, Inc . To facilitate the data 

preparation for the NAP program, the Minitrack preprocessor 

was modified to punch out the three control cards (Cotegories 20 I, 202, and 

999) required for egch data message. (For a description of the required 

control cards, see reference 4). 

Although nominally there is a single time associated with each 
pair of direction cosines, this is not so in practice. Differences 
arise due to different filter delays associated with each direction 
cosine and also due to different counter delays (reference 3). As 
the volume of printed output from the NAP program is proportional 
to the total number of different time points, the output volume 
could be reduced by one half, if each pair of direction cosines 
were adjusted to the same time point. This was done. The ad- 
justed time point was chosen as the average time of the two measure- 
ments. For a typical NAP run, this reduced the output volume by 
half. 

A modification was also made to output.on tape in a NAP input 
format, the fitted midpoint values of the original program. The 
idea behind this was that the fitted midpoint values could be used 
initially as an input to NAP, the "all data points” input being 
reserved for the final iterations through NAP. This in theory should 
have reduced the total computer time required. In practice, however, 
it was found that the "all data points" input when processed through 
NAP, provided the user with an extremely useful criterion for deter- 
mining the quality of the data in a Mini track message — the standard 
deviation of the error of all measurements within a message. For this 
reason, little use was made of the "fitted midpoint" input. 
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To help in the analysis of the results of the NAP program, a 
modification was made to print out a summary of the results of 
the preprocessing. This summary contained the following informa- 
tion for each station: For each data message the midpoint time, 
the three midpoint direction cosines and their rates of change, 
plus the prepass calibration constants. 

Finally, an error in the way the 2 Hz filter delay was being 
applied was corrected. 

A Program Error. Messages spanning midnight are not 
handled correctly by the program. A Minitrack message 
typically spaces a time interval of one minute or one 
half of a minute, so that this occurs very infrequently, 
but the program error should be corrected. 

The Minitrack preprocessor exists as a modified version of the original program 
developed by GSFC. Some of the information used in the original program is not 
used in its present form. The Input cards to the program, however, are the same 
as with the original program. Consequently, a large number of cards containing 
little information are required. 

The following information is required for operation of the program. 

Function Unit Description 

Card Input Data Set 5 Station Calibration Constants 

Card 1 Blank 

Card 2 Station Constants 

Card Column Word Definition Format 

1-6 STATIO Any 6 alpha- A6 

numeric 

characters 

7 Blank 


58 



Function 


Unit 


Description 


Card Column 

Word 

Definition Format 

00 

% 

«o 

KSTA 

Station number 

12 

10-13 

KFA 

Time adjustment 
EW-channel 

14 

, 14-17 

KFB 

Time adjustment 
NS-channel 

14 

18-20 


Blank 


21-24 

EWM 

Phase bias medium 
antenna, EW-channel 

F4.3 

25-27 

CLEWM 

Cable length in- 
equality, EW- 
medium channel 

F3.3 

28-31 

EWC 

Phase bias coarse 
antenna, EW-channel 

F4.3 

32-34 

CLEWC 

Cable length in- 
equality, EW- 
coarse channel 

F3.3 

35 


Blank 


36-39- 

EWFEQ 

Phase bias equatorial ~ F4.3 
fine antenna, EW-channel 

40-43 

EWFPO 

Phase bias polar' fi ne 
antenna, EW-channel 

F4.3 

44-46 


Blank 


47-50 

NSM 

Phase bias medium 
antenna, NS-channel 

F4.3 

51-53 

CLNSM 

Cable length in- 
equality NS-medium 
channel 

F3.3 

54-57 

NSC 

Phase bias coarse 
antenna, NS-channel 

F4.3 

58-60 

CLNSC 

Cable length in- 
equality NS-coarse 
channel 

F3.3 

61 


Blank 
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Function 


Unit 


Description 


Function 
Card Input 


Card Column 

Word 

Definition 

Format 

62-65 

NSFEQ 

Phase bias equatorial 
fine antenna, NS- 
channel 

F4.3 

66-69 

NSFPO 

Phase bias polar fine 
antenna, NS-channe 

F4.3 

1 

70-74 


Blank 


75-79 

DATE 

Date of calibration 

16 

Cards 3 thru 10 


Blank 


Cards 11-135 


15 sets of cards 2 thru 
10 are required. One 
set for each of the 
Minitrack stations 

Card 136 
Card Column 

Word 

Definition 

Format 

1-5 

KSAID 

Satellite ID 

15 

6-24 


Blank 


25-32 

FREQ 

Satellite transmitter 
frequency (MHZ) 

F8.3 

Card 137 


Blank card terminates 
data set . 



Unit Description 

Data Set 8 Data Selection Control Cards 

Card 1 80 alphanumeric characters used to 

write a message in the pre- 
processed printout 

Card 2 (more than one card can be used) 


Card Column 

Word 

Definition 

Format 

1-10 

now. 

Program will process 

215 


IHIGH 

messages between 




ILOW and IHIGH 



\ 
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Card Column Word 


Definition Format 


11- 15 IYEAR 2-digit integer for 15 

year of measurements 
e .g . 1966 = 66 

16-20 NEWARC Arc number assigned 15 

to measurements. If 
several ILOW, IH1GH 
cards define some arc, 

NEWARC should be 
blank on all but first 
card . 

21-25 NAPEND Normally left blank 

-1 = more data tape to follow 

0 = normal termination 

1 = no terminal record written 
(for tape addition) 

According to Ref. 4, the input 
tape to NAP should terminate 
with a negative number in the 
last record. However, with current- 
versions of NAP this is no longer 
necessary. If NAPEND is left 
blank, this terminal record is 
written on the tape. If on the lasty, 
"ILOW, IHIGH" card NAPEND 
is set equal to a positive non-zero 
integer, then the terminal record 
is not written. (The reason for this - 
option is to facilitate combination 
of two or more data tapes into a 
single tape). 

If the messages being processed 
are contained on more than one - - 
tape, then NAPEND should be 
set to a negative integer on the 
last "ILOW, IHIGH" card referring 
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Card Column Word 


Definition 


Format 


Function 
Tape Input 

Tape Output 

Tape Output 


to a tape that is to be followed 
by another tape . e .g . , if data 
messages to be processed are 
contained on three tapes such that 

FT09F001 defines tape 1, 
FT09F002 defines tape 2, 
FT09F003 defines tape 3, 

then the last "ILOW,IHIGH" 
card of tape 1 should have a 
negative NAPEND, the last 
"ILOW, IHIGH" card of tape 2 
should have a negative NAPEND,- 
the last "HOW, IHIGH" card of 
tape 3 should have either a blank 
or a positive NAPEND. 

Card 3, 4 Blank 

Each "ILOW, IHIGH" card should 
be followed by two blank cards. 

- - - The program actual ly has a restart -w 

capability. The blank cards must 
be inserted if the user does not 
wish to avail himself of the re- 
start facility. For use of the 
restart facility the user is referred 
to the program listing. Appendix 2.3 

Unit Description 

Data Set 9 Data tape(s) containing Minitrack messages to be 
processed , 

Data Set 12 Preprocessed output of "all data points" in a NAP-II 
input format 

Data Set 19 Preprocessed output of "fitted midpoint" in NAP-II 
input format 
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Function 

Unit 

Description 

Punched 

Data Set 7 

NAP-II control cards for each data message. 

Card Output 


These are .he start-stop times for that station 
data set (NAP-II cards group 201, 202, 999) 

Printed 

Data Set 6 

Raw Minitrack messages with each message fol- 

Output 


lowed by the computed direction cosines for each 
data point. (Direction cosines are also given in 
terms of Minitrack counts.) 


Data Set 1 1 

Printout of intermediate results obtained in pre- 



processing 


Data Set 13 

Input station constants (See data set 5) 

List of all messages that have not been preprocessed 
and reasons why. 

Pass summary information (number of messages per 

- 


station, time of last message for each station) 
Information required for restart. 


Data Set 14 

Summary of preprocessed results arranged in time 
sequence (station, time, direction cosines, 
direction cosine rates) 

- ' - 

Data Sets 

32 data sets — 1 per station.' Although there are 


20 thru 51 

only 16 stations the program treats a station operating 


in the equatorial mode as one station with station 
number KSTA (See data set 5) and the same station 
operating in the polar mode as a different station 
with station numbers KSTA+100. 

Each data set contains preprocessed results (for 
each station) arranged in time sequence (time, 
direction cosines, direction cosine rates, prepass 
calibration constants.) 

Sug gested Running Time 

For the 360/91 
For the 360/95 
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CPU time 11 minutes per 1000 messages 

CPU time 14 minutes per 1000 messages 



2.3 Optical Data 

2.3.1 Standard Optical Data Format 


The optical data used in this study was obtained from the NASA Space 
Sciences Date Center (NSSDC) at Goddard Space Flight Center. The data was 
on magnetic tape in the "GEOS format." A complete description of the data format 
is available from the NSSDC. The "GEOS format" for optical data is basically 
80 column card images stored on magnetic tape. Each card image contains the 
following. 


Card Column 
1-6 

7 

8 

9-11 

12-13 

14-18 

19-34 

35-53 

54-59 

60-71 

72-80 


Description 

Satellite identification 

Type of coordinates (RA and DEC) 

Observation identifier 

Timing standard deviation 

Time identifier 

Station number 

GMT of observation 

Observation data 

Date of plate reduction 

Code information as to processing 

Description of random error 


A separate data tape for each tracking network is normally obtained from the 
NSSDC. For example, in this study optical data from the STADAN network and 
the SAO network was to be used. Two data tapes were obtained, one for STADAN 
and one for SAO. In order to simplify the optical preprocessing, it was necessary 
to merge the two tapes. A special program was written to this task. The program 
listing is give., in Appendix A-2.4. The program merges two time sequenced optical 
data tapes into a single tape which will be the input data tape for the optical pre- 
processor. Data for each month is written on a separate file numbered sequentially 
on the same tape . 
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The following is required to execute the merge program: 

Function Unit Description 

Card Input Data Set 5 Tape and file number of output tape used for 


printout . 




Card 

Column 

Word 

Definition 

Format 

1-8 

TAPE 

8-character tape 
number 

A8 

9-10 


Blank 


11-15 

IFILE 

File number of first 
month of the 
combined tape 

15 


Tape Input Data Set 1 Data tape with STADAN stations. 

Tape Input Data Set 2 Data tape with SAO station. 

NOTE: The user should examine the printed output of the 

input data tapes to determine if there are dummy 
or non-essential data records on the tapes. 

Tape Output Data Set 3 Merged data output (Input tape to Optical 

Preprocessor^ JCL cards required for each 
month there is data. FT03F001 should have 
tape and file no. corresponding to card input, 
data set 5. The reason for writing a separate 
file for each month is to minimize tape search 
required in the optical preprocessor. 

Printed Data Set 6 Printed output of the merged output data tape 

Output (data set 3). Each page has a heading giving 

the year and month of the data and the tape 
and file number on which the data may be 
found on tape . 
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2.3.2 Optical Data Preprocessing 

The program assumes that the data on the input tape (data set 1) is arranged 
in time sequence. The program resequences the data such that the data for each 
station/orbit combination is sequential . An orbit is defined as consisting of TORBT 
seconds (45 minutes). Furthermore, each "station/orbit" is subdivided into passes 
(photographic plates) where each pass may be TPASS seconds (45 seconds) long. 

The data itself is converted to radians and output on data set 2 in the NAP input 
format. NAP control cards are punched (data set 7) for each pass. 

For SAO stations (station ID 29XXX), the input data is given in the 1950.0 
coordinate system. This is converted to the "true of date" coordinate system. The 
precession and nutation formulae used are based on those given in the Explanatory 
Supplement to the Astronomical Ephemeris and the American Ephemeris and Nautical 
Almanac (Her Majesty's Stationary Office, London, 1961). Also, the time of the 
SAO stations is given as Atomic time. This is converted to UTC using information 
provided by the user (data set 4). A further complication arises in the case of the 
SAO stations because, for those stations, the time associated with the data is the 
time of observation rather than the time that the signal left the satellite. An "r/c” 
correction should thus be applied to the timing of the orbit determination program. 
Unfortunately, NAP-II does not possess this capability. A temporary fix has, 
therefore, been made (subroutine SAOCOR) utilizing the fact that, nominally, 
the GEOS satellite always flashes on the even second. Since this is not exactly 
true, it is recommended that the user compare the output time for the SAO stations 
with those of other stations and make the appropriate adjustments on NAP "704" 
cards (see reference 4) . 

The following cards are required for running the optical preprocessor. 

Function Unit Description 

Card Input Data Set 5 Station and data set control cards. 
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Function' 


Description 


Uni t 

Station Cards: Up to 30 observing stations may be 

used. One card for each station. 

Card 


Column 

Word 

Definition 

Format 

1-5 

KSTA 

Station No. given on 
NAP 201/202 cards. 

15 

6-10 


Blank 


11-18 

STNAM 

8-alphanumeric character 
station name . 

A8 

19-20 


Blank 


21-25 

ISTA 

Station ID 

15 

26-30 

IEND 

Blank except for last 

15 


"station" card where a 
negative integer is used 
to indicate last station 
card . 

Data Set Cards: One card for each arc of data. 


Card 


Column 

Word 

Definition 

Format 

7-5 

NARC 

Arc number required for 
NAP control cards . 

15 

6-10 


Blank 


11-16 

IYMDB 

Year, month and day of 
beginning of arc 
(e.g. 660312). 

16 

17-19 

IHB 

Hour of beginning of arc . 

13 

20-22 

IMB 

Minute of beginning of 
arc . 

13 

23-32 

SECB 

Seconds of beginning of 
arc . 

D10.0 

33-38 

IYMDE 

Year, month and day of 

16 


end of arc . 
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Function 


Unit 


Description 


Card Input 


Tape Input 


Tape Output 

Printed 

Output 

Punched 

Output 


Card 

Column 

Word 

Definition 

Format 

39-41 

IHE 

Hour of end of arc . 

13 

42-44 

IME 

Minutes of end of arc . 

13 

45-54 

SECE 

Seconds of end of arc . 

D10.0 


Data Correction from A. I to UTC time. Required only 

Set 4 if SAO stations are to be processed. Any number 

of these cards may be used and the program will 
use a linear interpolation to any required time. 


At least two cards 

required . 


Card 

Column 

Word 

Definition 

Format 

1-2 

IYR 

Year 

12 

3-4 

IMO 

Month 

12 

5-6 

IDA 

Day 

12 

7-10 


Blank 


11-20 

COR 

A.I minus UTC on 
year, month, day 
(IYR, IMO, IDA) 

F10.5 


Data Set 1 


Data Set 2 
Data Set 6 
Data Set 3 


Data tapes to be processed (output from Merge or 
GEOS formatted data tape). If more than one 
data tape is used, then JCL will be modified to 
accommodate (e.g. FT01F001 for first tape, 
FT01F002 for second tape, etc.). 

Preprocessed data tape to be input to NAP-II. 

Print of preprocessed data. 

Printed output of input (data sets 4 and 5) . 


Data Set 12 
Data Set 7 


Summary of data . 

NAP Control Cards for timing/station information, 
NAP Cards 201, 202, 999 (see Reference 4); 


68 



2.4 REDUCTION OF MINITRACK DATA 
2.4.1 PRE NAP Program 

The processing of 20 days of Minitrack data requires about 5000 input cards. 
In order to reduce the number of cards that must be handled for a NAP-II run, the 
NAP-II control cards are written on magnetic tape. This "control card" tape is 
then updated and communicated to the NAP-II program via the PRENAP program. 
In this mode of operation, execution of the PRENAP program becomes the first 
step of a complete NAP-II run. 

The following is a listing of the program to copy the NAP-II control cards 
onto magnetic tape . 

Job Card 

// EXEC FORTRANG 
//SOURCE. SYSIN DD* 

IMPLICIT REAL * 8 (A-H,0-7) 

DIMENSION KAT(2),KEY(10),DATA(2) 

100 FORMAT (13, 12,A8, 101 3, D22 .8, D 15 . 8) 

DO 300 1=1, 10000 

READ (5, 100,END=400) KAT,XLABEL, KEY, DATA 
WRITE (9) KAT, XLABE L, KEY, DATA, I 
300 CONTINUE 
400 ENDFILE 9 
REWIND 9 

200 FORMAT (1X,3I8,A8, 1013, 2D19.8) 

DO 330 J=1 , 10000 

READ(9,END=500) KAT, XLABEL, KEY, DATA, I 
WRITE (6, 200) I, KAT, XLABE L, KEY, DATA 
330 CONTINUE 
500 STOP 
END 
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// EXEC LOADER, PARM='MAP, CALL, SIZE=100K', REGION .GO=l 10K 
//GO.FT09=001 DD UNIT=9TRACK,LABEL=(4,BLP),DSN=MRGDAP,DISP=(NEW,KEEP) 
// DCB=(RECFM=VBS / LRECL=80 / BLKSIZE=3204) VOL=SER=34517C 
//GO.DATA5 DD* 

Once the NAP-II cards are written on tape, the function of the PRENAP 
program is to update the "cards" on tape and to convert the binary "cards" tape to 
a format acceptable to the NAP-II program. The PRENAP program is listed in 
Appendix A-2.6. At the end of the program listing is an example of the JCL cards 
needed and a sample set of update cards. 

The program is normally loaded as a binary object deck (Appendix A-2.6 
JCL cards). The function of the peripheral equipment needed to execute the 


PRENAP program are: 



Function 

Unit 

Description 

Input Tape 

Data Set 9 

NAP-II control cards on magnetic tape. 
(Output from program given above).. 

Disk File 

Data Set 10 

Used as an intermediate scratch tape. 

Disk File 

Data Set 12 

Used to pass "updated" NAP-II control cards 
to NAP-II program for execution. >_■ : 

Card Input 

Data Set 5 

Updater cards for NAP-II control cards. 

If no update is required, data set 5 should 
be left empty. There would probably be 


no update the first time NAP-II is run for 
a particular job. The program prints out 
instructions on how to update "cards". 

2.4.2 Use of NAP-II Program 

The NAt'-II program is a very flexible analysis tool providing the user with 
considerable freedom for designing and executing data reductions. Consequently, 
the user must have a good understanding of the program uses to attain efficient 



utilization . This can only be done through use of the program and its various options. 
The functions of the program are discussed in reference 4. 

The type of information required to execute NAP-II is given as follows.- 

Type of Information Categories 

General and planetary information 101, 102, 103, 151 

Comments 150 

Station survey 301, 302 , 303 

Totally stable parameters 601, 602 


Measurement definitions 

General arc information 
Arc stable parameters 
General pass information 
Pass stable parameters 
Pass comments 


701, 702, 703, 704 
(Sets 0 only) 

104, 205, 206 

601, 602 

201, 202 , 203 , 204 
601, 602 
152 


Overrides for measurement definitions 701, 702, 703, 704 

for a pass (Sets 1 only) 

End of pass, arc, or all control data 999 

The station numbers (category 300) are preassigned by the Minitrack and the op- 
tical data preprocessors, depending on the order of the Minitrack Station Calibration 
Cards and the Optical Stations Control Cards. This controls both the station ID num- 
bers and the user assigned station numbers. These two programs also provide the timing 
information cards (category 200 and 999) as part of their output. The other cards must 
be set up by the user. The information required by the various groups of cards is given 
in reference 4. 


Once a "master" set of cards has been set-up, subsequent runs can be made with 
new data sets by changing the following quantities, 

© State vectors - Category 601, 205 

© Greenwich Hour Angle - Category 206 
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o Lunar and Solar Coordinates - Category 104 
© Timing information - Output from data preprocessors. 

These cards can be changed via the PRENAP program. 

To run the NAP program, it is essential to have a reasonably good estimate of the 
state vector at the initial epoch. Before the user tries to solve for any parameters, he 
should make an initial run, which will allow the POSTNAP program to edit out ob- 
viously bad data and correct the lobe assignments made by the Minitrack preprocessor. 
Once this has been done, the user can go ahead and solve for any parameters he wishes. 
As the user improves his estimate of the initial state vector, he may wish to modify his 
definition as to what constitutes bad data. 

The NAP-II program outputs the following information on disk or tape: 

measurement discrepancies (FT37F001) 

current estimates of error model parameters (FT37F002) 

data tape with the mid-data point for each 
pass (FT37F003) . 

For multi -arc data processing, it is recommended that the user switch from using the 
data tape output from the Mini track preprocessor to the "single data point per pass" 
data tape. When using the single data point prepass, the Category 704 cards have to 
be deleted, because the biases have already been applied. The run time advantage of 
using one data point instead of all data points (which may be as many as 155) is con- 
siderable. The loss of accuracy is negligible. 

The recommended procedure for reducing Minitrack data is summarized in 
Section 2.6. 

There were changes made to the NAP-II program to effect the procedure for re- 
ducing Minitrack data. The need for some of these changes was discussed in the preface 
of this report. Without going into great detail, the changes to the particular subroutines 
are as follows: 
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© Changes to Subroutine ENEXPS 

This subroutine computes two sets of series: (1) the power series ex- 
pansions of n-heavenly bodies, and (2) the power series expansion 
defining whether or not the satellite is exposed to solar radiation. 

The subroutine was completely reprogrammed, making the new version 
seven times faster than the old. One of the factors contributing to 
the increased computational speed of the new program is that it 
computes the n(n-l)/2 different distances between n bodies, and no 
more. The old program somehow managed to compute more distances 
than this so that some distances had to be identical . 

The old program had a programming error in the computation of the 
"shadow" series so that the two programs do not agree as far as 
shadow series expansion is concerned. The FORTRAN listing of the 
subroutine is given in Appendix A-3. 1 . 

© Changes to Subroutine EXPAND 

Subroutine EXPAND computes two sets of series: (1) the power series 
expansion of the acceleration due to the full gravity field, ex- 
pressed as a sum of spherical harmonics, of the main attacking body, 
and (2) the power series expansion of the acceleration due to atmos- 
pheric drag. The part dealing with the gravity field was completely 
rewritten. The mathematical analysis for this is given in reference 5. 
The FORTRAN listing of the subroutine is given in Appendix A-3. 2. 

© Changes to Subroutine VARIEQ 

Subroutine VARIEQ computes the two sets of series: (1) the power 
series expansion of the variational matrix due to the full gravity field 
of the main attracting body, and (2) the power series expansion of 
the variational matrix due to atmospheric drag. The FORTRAN 
listing of the subroutine is given in Appendix A-3. 3. 
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© Changes to Subroutine FINALP 

A small change was made in this subroutine to output on 
disk (data set 10) the computed values for the error model 
terms. Note that the call to this routine precedes the final 
call to subroutine RESID so that the quantities output by this 
routine on data set 10 precede those output by RESID. The 
FORTRAN listing is given in Appendix A-3.4. 

© Special Subroutine PRTIAL 

The standard subroutine PRTIAL computes the discrepancies 
between observed and predicted measurements and also the 
partial derivatives of the predicted measurements with respect 
to those error model terms that affect the predicted measure- 
ments. The computed values are then stored on a disk file 
(ISFILE). Later in the program subroutine SOLVER forms the 
normal equations from this disk file and solves the equations 
to give the computed values for the error model terms. The 
quantities on the disk file are also used by subroutine RESID 
which prints the measurement residuals. (The residuals are up — 
dated discrepancies between observed and predicted measure- 
ments, where the predicted measurements reflect the latest 
values (from SOLVER) of computed error model terms. The 
quantities involved in this computation are the discrepancies 
and partial derivatives from subroutine PRTIAL, which are 
stored on disk file, and the difference between the computed 
values for the error model terms and the values used by sub- 
routine PRTIAL for the predicted measurements). 

Obviously, the more data points that are stored on ISFILE the 
longer it will be. For a one-week long arc of Minitrack data, 
ISFILE is about 3,000,000 bytes long (20 cylinders of disk space) 
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To print the final residuals would require a further 9,000,000 
bytes (60 cylinders of disk space). Together with the remaining 
disk requirements of the NAP program this puts a severe strain 
on the total 360 system. 

To make it practicable to use the NAP program for multi-arc 
cases, several changes were made to subroutine PRTIAL. 

Changes to Subroutine PRTIAL. For each measurement all 
discrepancies for a pass (Minitrack message) are averaged 
to give the mean measurement discrepancy. The standard 
deviations are also computed and any discrepancies differing 
from the mean by more than three standard deviations are 
rejected. (The process is iterative, since as soon as a dis- 
crepancy is rejected a new mean and standard must be 
computed). The mean discrepancies are then output on 
ISFILE together with the appropriate partial derivatives. 

The time associated with the mean discrepancies is also 
computed. These measurements are written on the same 
tape on disk file as the discrepancies (FT37F001). They are 
also written on disk (FT35F001) in a NAP observation tape 
format. The standard deviations and the number of measure- 
ment points (not rejected) are written on disk (FT36F001). 

The FORTRAN listing of the subroutine is given in Appendix 
A-3.6. 

© Changes to Subroutine RESID 

There-are two subroutines RESID that can be used in NAP-II. These 
are identified as STANDARD Subroutine RESID and SPECIAL Sub- 
routine RESID. The Special subroutine was written to make multi- 
arc runs practical in terms of run time. 
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Standard Subroutine RESID. Three sets of changes were 
made to this subroutine. The new subroutine will work 
only if there are two or fewer kinds of measurements per 
station-pass. The first change involved writing the printed 
(as on data set 6) pass summary on disk (data set 10) for 
use by some postNAP program, and also on data set 34, 
which is a printed pass summary. The second change in- 
volved plotting the measurement residuals in addition to 
giving the numerical values. The third change involved 
the computation of the correlation coefficient between the 
errors of the two measurements. The reason for computing 
the correlation coefficient was that when the plotting 
facility was added to this subroutine, it looked as if 
there was some correlation between the two measurement 
errors. However, in most cases the correlation coefficient 
is very small . The FORTRAN listing of the subroutine is 
given in Appendix A-3.5 . 1 . 

Special Subroutine RESID. The computed values of the 
error model terms are copied from disk (FT1 0F001) and 
written on tape (FT37F002). The measurement residuals 
are computed as before. However, since there now is only 
one measurement time point per pass, the standard deviation 
cannot be computed. Instead, the standard deviation computed 
in subroutine PRTIAL is copied from disk (FT36F001) as is the 
number of measurement points. This information is written on 
the same tape as the new values for the error model terms. 
Finally, the new observation tape is copied from disk (FT35F001) 
to tape (FT37F003) . The FORTRAN listing is given in Appendix 
A-3.5. 2. 
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© Changes to Subroutine MESOLD 

Terms were added to Subroutine MESOLD to reflect the complete 
Minitrack geometrical error model. The following should be added 
to the NAP-II User's Guide (Reference 4) to be consistent with these 
changes (Appendix IV, Page IV-2). 



where measurement code 6 corresponds to direction cosine alpha 
(or t), and measurement code 7 corresponds to direction cosine beta 
(or m), and n = yf\~ (-L 2 + m 2 ) 

The 12th error model term (TM) is not part of the geometrical error 
model but corresponds to a timing bias. 

RE and RP are rotation terms about the local vertical of the equatorial 

and polar m measurements, respectively. RL and RM are rotations or tilt about 

the t or m measurements. 

2.4.3 The NAP-II Minitrack Post- Processing Program 

The post- processing program summarizes the results of a NAP-II Minitrack run and 
prepares the NAP-II "input cards" for the next iteration. Any lobe assignment errors of 
the Minitrack preprocessor can also be corrected through the NAP-II "input cards". 

In order to execute the post-processing program, the NAP-II program must be a 
version employing the special subroutine PRTIAL and RESID given in Appendix A-3.6 
and 3.5.2, respectively. 

The program assumes that the station numbers assigned by NAP are odd for equatorial 
stations and even for polar stations and that measurement numbers are odd for "L" measure- 
ments and even for "M" measurements. (See Appendix A-4.) 
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The NAP-II control cards must be in a particular order required by the POSTNAP 
program. In general, error model terms should be setup in order of totally stable, arc 
stable and pass stable. The error model parameter numbers (Category 601, Key 7), must 
be in ascending and consecutive order. Cards must not be repeated. 

This is not- a restriction on the NAP-II program but a restriction imposed by POSTNAP. 


Function Unit Description 

Card Input Data Set 5 Only one card required. 

Card 

Column Word Definition 

1-5 IOPT Determines which NAP "701" (con- 

tinuation 1) cards are to be output. 

If measurements were previously 
edited out of the solution and no 
value for the pass standard deviation 
is available, the measurements are 
edited out again.- IOPT can assume 
values of 0, -1 and +1 . In all three 
cases the above applies . In all cases 
key 7 of 701 card is set to "1 • 

= 0, measurements are edited out if mean 
pass measurement error exceeds MAXERR 
or if pass standard deviation exceeds 

... . SDMAX. If measurements are edited 

- - --- - 2 • out, DATA 2 of 701 card is set to 1 .D15. 

= -1, measurements are edited out if mean 
pass measurement exceeds MAXERR or 
if no pass standard deviation is avail- 
able or if pass standard deviation exceeds 
SDMAX. If measurements are edited, 
DATA 2 is set to 1 . D 1 5 . 

— ... = +1, same as -1 case except if measure- 

ments are not edited out, they are 
assigned "a priori" sigma of 10 times 
pass standard deviation. If measurements 
are edited out, DATA 2 is set to T.D7 
times the pass standard deviation . 
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Function 


Unit 


Description 


Card 


Column 

Word 

Definition 

6-10 

MAXERR 

Maximum mean measurement error 
in Minitrack counts. 

11-20 

SDMAX 

Maximum pass standard deviation 
(unsealed) . 

21-30 

FREQCY 

Satellite transmitter frequency in 
MHz. 


Note on the avai labi lit y of a Pass Standard Deviation . 

A pass standard deviation can only be computed if there is more than one measure- 
ment per pass. If IOPT = 1, the pass standard deviation is saved on a "701" card. 
(If key 7 =0, it is equal to DATA 2 times 1 . D— 1 ; if Key 7 = 1, it is equal to 
DATA 2 times l.D-7). Use is made of this facility when the data processing 
switches from using all measurements for a pass to the single new measurement 
computed in special subroutine PRTIAL. 

N ote on data editing . 

A bad pass is defined as one whose mean error exceeds MAXERR or whose 
standard deviation exceeds SDMAX. Bad passes are edited out of the solution 
in all cases. If no standard deviation is available, then in the case of Iopt = -1 
the pass will be edited out of the solution. (It is recommended that Iopt = -1 
only be used if there are more than one data point per pass. The non-availability 
of a pass standard deviation may in that case be the result of the pass not meeting 
the elevation requirements. NAP in that case would output a zero mean pass error) 

Iopt = 1 should be used if the user wishes to weight the measurements of a pass in 
inverse proportion to their standard deviation, in other words, a noisy pass would 
be given less weight than one that is less noisy. 


79 




Function 


Unit 


Description 


Disk Input 
(on tape) 

Data Set 27 

Disk inp d 

Data Set 10 

(on tape) 

File 1 
(FT10F001) 

Disk input 

Data Set 10 

(on tape) 

File 2 
(FT10F002) 

Disk output 

Data Set 10 

(on tape) 

File 3 
.(FT10F003) 

Printed output 

Data Set 6 

Printed output 

Data Set 20 


Printed Output Data Set 31 


NAP input cards used in making NAP run. This 
would normally be on a disk file created by the 
PRENAP card updating program . 

Output from NAP data set 37 (FT37F002) . The 
computed values of error model parameters. 

Output from NAP data set 37 (FT37F001). Current 
measurement discrepancies. 

Output of sequenced NAP "input cards" for use 
by the PRENAP card updater program on the next 
iteration . 

Prints the following 

© Some instructions on the use of the program, 
e Pri ntout of input cards (data set 5) . 

© Printout of NAP "output cards" (FT10F003) 
for the next iteration. The first column of 
this list are the sequence numbers . - - 

The following is printed out for each station/pass. 

© Time of midpoint of pass . 

© Mean measurement error (This is the mean 
discrepancy computed in special subroutine 
PRTIAL and equals the mean residual error 
computed in subroutine RESID when the 
solution converges). 

© Standard deviation of the measurements for a pass. 
For only one measurement the standard deviation . 
is equal to zero. . - 

© For each measurement, the time relative to mid- 
point and the error relative to mean pass error. 
This error is also plotted. 

Printout of the bias values used in the Mini track 
— preprocessor (built into the program) and the new 
bias. New biases are obtained by converting the 
bias obtained by NAP to Minitrack counts and 
adding the result to the preprocessor bias. 
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Function 


Unit 


Description 


Printed output 


Data Set 31-60 A summary for each station is written on a 
different data set. Each pass contains 

© Time of midpoint measurement 
© Arc number 
© Pass number 

© Mean error in Minitrack counts 
© Whether ( ) or not (#) the measurement 
met elevation requirement given on 
NAP "702" card 

© Whether ( ) or not (*) measurement con- 
tributed to present solution 
o Measurement standard deviation in Minitrack 
counts 

e Bias correction made to the measurement for 
the next iteration expressed in Minitrack 
lobes (in 1000 Minitrack counts) . The 
errors being corrected are due to wrong 
lobe assignments in the Minitrack preprocessor 
and are corrected via NAP "704" cards 
e Number of measurement points per pass on the 
data tape 

o Number of measurement points deleted by 
special subroutine PRTIAL 
© The mean error rate 

- © The reduction in standard deviation,- expressed 

as a percentage, that would result if the 
standard deviation were computed about 
— - the line defined by the mean error and 

mean error rate and not the mean error alone 
© The computed measurement. The predicted 
(NAP) measurement plus the mean pass 

measurement error, (computed in special 

subroutine PRTIAL). 

e Whether ( ) or not ■(*) the measurement will 
- - ... contribute to the solution of the next 
iteration. 
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2.5 RESULTS OF MINITRACK DATA REDUCTIONS USING NAP-II 

2.5.1 Estimation of Calibration Parameters for the Fort Myers Minitrack Station 
Using the GEOS-I Satellite and Optically Determined Reference Arcs 

A number of GEOS-A arcs determined from optical observations were used as 
reference arcs. These arcs had previously been obtained by DBA Systems, Inc. for 
work on another contract. Since optical observations are considerably more accurate 
than Minitrack, the optical reference arcs could be considered perfect. The NAP 
program was then used to recover measurement biases and measurement scale factor 
errors. The Minitrack measurements for each pass were weighted in inverse pro- 
portion to the variance of the measurement error for each pass (Minitrack message) . 

The results are shown in Tables 2 and 4. Tables 1 and 3 give the actual 
measurements. Table 5 is a comparison of calibration results. Table 6 tabulates 
some observed frequencies. It shows that the fluctuation in the GEOS-A transmitter 
frequency is very small . Table 7 shows the effect of including the Doppler effect 
and wave propagation time delay in the computations. The total effect can be seen 
to be of the order of 1 Mini track count. 
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TABLE 1 


POLAR ARRAY. Observed values of direction cosines, direction cosines 
rates and prepass internal bias (KS2) . 
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TABLE 2 


POLAR ARRAY. Measurement Errors. 


. , Standard Deviation Of 


Mean Measurement Error For Pass 


^ Measurement Error About Mean 
i , For Pass 


a! 

am 

A1 

Am 

.000072 

.000147 

-.000001 

-.000135 

.000036 

.000034 

-.000038 

-.000097 

.000749 

.000660 

-.000164 

-.000034 

.000068 

.000072 

- .000027 

.000085 

.000040. 

.000063 

.000009 

.000033 

.000044 

.000056 

.000116 

-.000005 

.000239 

.000141 

-.000042 ' 

-.000023 

.000261 

.000353 

.000033 

-.000057 

.000044 

.000032 

.000155 

-.000068 

.000071 

.000081 

.000048 

.000125 

.000133 

.000173 

.000112 

.000029 

.000042 

.000048 

.000011 

.000005 

.000035 

.000034 

.000003 

.000098 

.000044 

.000057 

-.000057 

.000136 


.000038 

* 

.000036 

.000052 

.000048 

.000087 

.000047 

.000085 

.000092 

.000130 

.000118 

.000029 

.000032 

.000003 

.000018 

.000039 

.000034 

- .000046 

-.000005 

.000021 

.000028 

-.000062 

-.000066 


MINITRACK 
COUNTS 
A! Am 



EW - BIAS: 137 MINICTS - .000021 = 136 MINICTS 
NS- BIAS: 431 MINICTS +.000071 = 435 MINICTS 
EW- SCALE FACTOR: -.00022 
(Base Line is Shorter Than Assumed.) 
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ARC Number 


TABLE 3 


EQUATORIAL ARRAY. Observed values of direction cosines, 
direction cosine rates and prepass internal bias (KS2) . 



Direction Cosines 

h- 

2 

O 

DAY HRS. MIN. I m n 


Direction Cosine Rates 


| m n 


K 

o 

o 

I 


CN 

to 


o 

o 

LO 

I 


J , 

CN 
t n 
'iZ 




-.04 -.65 .76 

-. 03 - -.90 .43 


.00 .85 .53 
.00 .57 .82 
.01 _ .21 .98 
.00 .73 .68 
.00 .52 .85 


.0020 

.0022 

.0028 

.0023 

.0016 

.0019 

.0014 

.0014 

.0013 

.0014 


-.0007 .0013 

-.0013 .0017 

-.0042 .0011 

-.0027 -.0022 
-.0009 -.0018 2 

-.0008 .0013 1 -1 

.0014 -.0010 6 -42 
.0022 -.0005 5 -51 

.0009 - r .0010 16 14 

.0016 --.0010 18 i 15 
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TABLE 4 


EQUATORIAL ARRAY. Measurement Errors. 

1 — 

Standard Deviation Of Mean Measurement Error For Pass 


Standard Deviation Of 
Measurement Error About 
Mean For Pass 


MINITRACK 

COUNTS 


1 

.000191 

.000213 

000 1 5 1 

- 000080 

7 

-4 





2 

.000081 

.000053 

-.000013 

-.000050 

-1 

-2 

4 

.000033 

.000040 

.000062 

.000025 

3 

. 1 

13 

. .000189 


.000157 


7 

* 

15 

.000657 

.000757 

.000192 

-.000397 

9 

-18 

17 

.000083 

.000076 

.000056 

.000006 

3 

0 

20 


.000041 

* 

-.000007 

* 

0 

21 

.000031 

.000041 

-.000054 

-.000040 

-2 

-2 

23 

.000071 

.000043 

.000113 

.000012 

5 

1 

25 

; 

.000034 

.000039 

- .000045 

.000038 

-2 

2 


EW - BIAS: 957 MI NIC TS +.000026 = 958 MINICTS 
NS - BIAS: 988 MINICTS - .000544 = 988 MINICTS 

NS - SCALE FACTOR: -.00049 
(Base line is shorter than assumed.) 
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TABLE 5 


COMPARISON OF CALIBRATION RESULTS 



* The DBA computed calibration values are the ones shown in Tables I and II. 
They are based on observations from 4 January to 2 April 1966. 

The RCA results have been verbally communicated by Mr. Jerry Casto and 
are based ^n satellite observations covering approximately the same time 
period as those used by DBA. 
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TABLE 6 


OBSERVED GEOS - A FREQUENCIES 
(Source: Harry Pritchard, GSFC) 


Scale Error 
Relative To 
Nominal Value 


Date (1966) 

Value (MHz) 

(136.83 MHz) 

1 January 

136.830360 

26 x 10* 7 

5 January 

136.830350 

26 x 10' 7 

3 February 

136.830038 

3 x 10~ 7 

8 February 

136.829 675 

-24 x 10~ 7 

17 March 

136.830540 

I 

o 

X 

o 
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TABLE 7 


CHANGE IN PREDICTED MEASUREMENTS DUE TO INCLUSION OF 
THE DOPPLER EFFECT AND WAVE PROPAGATION TIME DELAY 



POLAR ARRAY 


• EQUATORIAL 

ARRAY 

ARC No. 

10 s A 1 

10 s Am 

ARC No. 

10 s A I 

10 s Am 

1 

-14 

19 

1 

-15 

17 

2 

-14 

20 

2 

-15 

18 

3 

-10 

21 

4 

-13 

20 

5 

-12 

21 

13 

-11 

X 

7 

-11 

21 

15 

-11 

21 

9 

-11 

21 

17 

-14 

18 

10 

-14 

18 

20 

X 

-17 

11 

-10 

21 

21 

-11 

-18 

12 

-12 

20 

23 

-13 

-16 

15 

-11 

21 

25 

-12 

-17 

16 

-14 

18 



• 

17 

-14 

19 




18 

-13 

20 




19 

-11 

20 




20 

X 

-18 




21 

-11 

-18 




22 

- 8 

-20 




23 

-12 

-17 




24 

-11 

-18 




25 

-11 

-17 





POLAR ARRAY: I MINITRACK Count = .0000175 

EQUATORIAL ARRAY: 1 MINITRACK Count = .0000216 
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2.5.2 The Calibration of all Minitrack Stations Using GEOS-A Minitrack Measurements 

The procedures followed in processing the Minitrack data are described in Section 
2.6 of this report. The four arcs described in Table 15 were at first processed separately 
A new data tape was created for each of the four arcs (see special subroutines PRTIAL 
and RES1D). The arcs were then combined in pairs to give two, two-arc runs, and 
finally all the data were combined into a single four-arc run. Three sets of runs were 
made under slightly varying conditions (see Table 15): In the first run scale factors 
and antenna rotation terms were recovered in addition to measurement biases and the 
four state vectors. The second run was similar to the first except that no rotation terms 
were considered. This resulted in very slight differences in the solutions. The third run 
was similar to the first but in this run ionospheric corrections were made using the 
ionospheric model developed by DBA under Contract NAS511730. In all runs, data 
were edited out of the solution using the criteria listed in Tables 12 and 13. It should 
be mentioned that when the runs were made, no ionospheric data could be obtained for 
1966. The 1971 val ues were therefore substituted. This may account for the relatively 
bad results obtained for run 3. Tables 8 through 14 refer to run 1 . Table 8 is a com- 
parison of calibration results. Table 9 lists the scale factor errors and antenna align- 
ment errors. It can be seen that the scale errors are predominantly negative. This is 
the kind of effect that would be expected if no ionospheric corrections were made, i.e., 
it seems likely that the scale factors have absorbed some of the effect of ionospheric 
refraction. When ionospheric corrections were made (run 3), the recovered scale fac- 
tors were indeed reduced in magnitude (the actual numbers are not included in this 
report. 

Table 10 is a listing of the number of observations for each station. Table 11 
gives the RMS values of the observed residuals. It is interesting to note that the Polar 
Array, which has a larger baseline, shows a larger RMS error than the Equatorial Array 
when the errors are expressed in Minitrack counts, but more or less the same RMS error 
in terms of angular error. From this, it could be inferred that the errors are caused less 
by the electronics of the observing stations, which measures Minitrack counts, than 
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by general uncertainties in conditions outside the observing station. This, in general, 
is probably true. However, it has been observed that large variations in the pre- 
calibration biases, which is part of the Minitrack message, have been accompanied 
by large residuals in the orbital fit. The Woomera Polar Station was particularly 
bad in this report. The measurement residuals and precalibration constants for that 
station are listed in Table 14 for illustrative purposes. Tables 12 and 13 are self- 
explanatory, as are Tables 15.1 through 15.4. It should be emphasized that the 
results given in Tables 15 are the most important results of this report. 

In all orbital computations the Earth model used was that given in Appendix A-l 
The affect of solar radiation pressure was not taken into account, and neither was 
the effect of the Earth's precession and nutation. 

In all three runs the measurements were weighted proportionately to the standard 
deviation of the measurements within each message. 
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TABLE 8 

Comparison of Calibration Results for 


Equatorial and Polar Station Biases 


Station 

This Progra 

m * 

Aircraft 

65.12.31 

Aircraft 

66.02.25 



EW 

NS 

EW 

NS 

EW 

NS 

Fort Myers (E) 

956 

981 

954 

988 

953 

984 

Fort Myers (P) 

135 

435 

138 

435 

137 

435 

Quito (E) 

817 

036 

812 

036 

837 

030 

Quito (P) 

756 

000 

760 

002 

763 

006 

Lima (E) 

949 

252 

955 

263 

955 

259 

Lima (P) 

023 

881 

028 

882 

028 

882 

Santiago (E) 

078 

057 

073 

060 

076 

059 

Santiago (P) 

950 

972 

943 

967 

943 

963 

New 

Foundland (E) 

930 

937 

934 

935 

934 

935 

New 

Foundland (P) 

170 

920 

171 

919 

171 

919 

Winkfield (E) 

071 

973 

067 

976 

067 

976 

Winkfield (P) 

868 

001 

870 

011 

870 

011 

Johannesburg (E) 

048 

937 

054 

932 

054 

932 

Johannesburg (P) 

084 

876 

. 090 

871 

090 

871 

Blossom Point (E) 

948 

094 

948 

096 

948 

096 

Blossom Point (P) 

081 

887 

082 

889 

082 

889 

College (E) 

955 

807 

961 

804 

960 

804 

College (P) 

846 

932 

854 

916 

859 

919 

Mojave (E) 

030 

576 

017 

579 

029 

575 

Mojave (P) 

079 

123 

062 

129 

079 

121 

Grand Forks (E) 

068 

951 

060 

954 

061 

951 

Grand Forks (P) 

602 

001 

598 

990 

598 

990 

Woomera (E) 

020 

944 

031 

938 

031 

938 

Woomera (P) 

934 

118 

948 

~ i}8 

948 

118 

* Time span covei 

red was 66. 1 

.1 - 66.2.3 

and 66.3. 

12 - 66.3. 

19. 
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Table 9 


Scale Factor Errors and Ante nna Alignment Errors In Minitrack Counts 

1000 Equatorial Minitrack Counts = .0216 radians 
1000 Polar Minitrack Counts = .0174 radians 

Equatorial Array Polar Array 

Scale Error Antenna Rotation Scale Error Antenna Rotation 


Fort Myers 

-14 

4 

-15 

2 

Quito 

9 

-4 

- 6 

-2 

Lima 

-24 

-4 

30 

-3 

Santiago 

- 2 

-8 

- 8 

-2 

New Foundland 

- 8 

0 

-13 

-2 

Winkfield 

- 4 

-3 

- 8 

5 

Johannesburg 

- 6 

-1 

-14 

-4 

Blossom Point 

-11 

6 

- 1 

0 

College 

-10 

8 

- 4 

0 

Mojave 

- 3 

-3 

- 7 

4 

Grand Forks 

8 

-5 

- 8 

-9 

Woomera 

-11 

-3 

- 8 

-4 

Scale Error 

The predicted measurement is 
equal to the "scale error" mu 
measurement. 

in error by an amount 
Itiplied by the actual 

Antenna Rotation 

The predicted measurement is in error because the 
antenna baselines are rotated anti-clockwise by the 
"antenna rotation". 


In the equatorial mode the scale error only affects the m measurement and the antenna 
rotation only the I measurement. 

In the polar mode the scale error only affects the I measurement and the antenna 
rotation only the m measurement. 
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Table 10 

The Number of Minitrack Messages Used in the Orbit Determination 


and the Determination of Biases, Scale Factors and Antenna Rotations 


Equatorial 

Array 

Direction Cosine 

1 

Total 

Direction 

Cosine 

m 

Total 

Arc 1 

Arc 2 

Arc 3 

Arc 4 

Arc 1 

Arc 2 

Arc 3 

Arc 4 

Fort Myers 

4 

4 

11 

12 

31 

4 

4 

11 

12 

31 

Quito 

1 

- 

2 

4 

7 

1 

- 

- 

3 

4 

Lima 

- 

3 

- 

2 

5 

- 

3 

2 

- 

5 

Santiago 

1 

1 

3 

3 

8 

1 

1 

3 

3 

8 

New Foundland 

17 

19 

26 

31 

93 

16 

19 

26 

31 

92 

Winkfield 

12 

- 

9 

11 

32 

12 

- 

15 

11 

38 

Johannesburg 

- 

1 

3 

2 

6 

- 

1 

3 

2 

6 

Blossom Point 

12 

9 

20 

16 

57 

12 

9 

21 

16 

58 

College 

12 

12 

18 

17 

59 

13 

12 

17 

17 

59 

Mojave 

12 

11 

17 

14 

54 

12 

11 

17 

14 

54 

Grand Forks 

11 

14 

22 

7 

54 

12 

14 

22 

7' 

55 

Woomera 

1 


1 

3 

5 

— 

- 

2 

3 

5 

Total 

83 

74 

132 

122 

411 

83 

74 

139 

119 

415 
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Table 10 (Continued) 

The Number of Minitrack Messages Used in the Orbit Determination 
and the Determination of Biases, Scale Factors and Antenna Rotations 


Polar 

Direction 

Cosine 

1 


Direction 

Cosine 

m 


Array 

Arc 1 

Arc 2 

Arc 3 

Arc 4 

Total 

Arc 1 

Arc 2 

Arc 3 

Arc 4 

Total 

Fort Myers 

11 

9 

14 

19 

53 

10 

9 

14 

19 

52 

Quito 

4 

3 

5 

5 

17 

4 

4 

5 

5 

18 

Lima 

5 

5 

5 

4 

19 

5 

6 

6 

6 

21 

Santiago 

4 

2 

5 

5 

16 

4 

3 

5 

4 

16 

New Foundland 

9 

14 

14 

17 

54 

9 

14 

14 

18 

55 

Winkfield 

3 

- 

2 

3 

8 

3 

- 

2 

2 

7 

Johannesburg 

2 

3 

7 

7 

19 

2 

3 

7 

8 

20 

Blossom Point 

13 

9 

15 

20 

57 

12 

9 

15 

20 

58 

College 

7 

5 

9 

8 

29 

7 

6 

9 

7 

29 

Mojave 

14 

11 

13 

20 

58 

14 

11 

13 

20 

58 

Grand Forks 

4 

7 

18 

10 

39 

4 

7 

18 

10. 

39 

Woomera 

4 

2 

3 

7 

16 

_3_ 

3 

5 

7 

18 

Total 

80 

71 

110 

125 

385 

77 

75 

112 

124 

389 
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Table 1 1 

RMS Value of Residuals^in Minitrack Counts 



Equatorial 

Array 

Polar Array 


Station 

1 

m 

1 

m 

Fort Myers 

2.7 

5.6 

5.1 

4.5 

Quito 

3.2 

8.9 

7.4 

5.2 

Lima 

3.7 

9.3 

6.6 

5.7 

Santiago 

4.7 

3.8 

7.2 

7.0 

New Foundland 

4.9 

3.3 

7.1 

5.3 

Winkfield 

7.3 

3.5 

4.5 

5.4 

Johannesburg 

3.6 

2.8 

4.5 

4.0 

Blossom Point 

4.8 

3.4 

4.8 

3.2 

College 

3.7 

5.9 

5.8 

8.4 

Mojave 

3.1 

5.3 

5.2 

5.1 

Grand Forks 

4.9 

4.2 

6.0 

7.1 

Woomera 

4.1 

11.8 

7.2 

8.1 

All ** 

4.3 

4.7 

5.9 

5.6 

*The residuals are 

the differences between 

observed and predicted values 

The predicted val 

ues are based on 

i the computed orbit. 



*#It is of interest to note that the four RMS residuals for all stations expressed as 
unsealed quantities assume the values: 

.000093, .000102, .000102, .00097. 

The same residuals expressed as equivalent angular discrepancies at zenith 
assume the values: 19", 21", 21", 20", i.e., values around the quoted 
Minitrack accuracy of 20". 
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Table 12 

Total Number of Minitrack Messages 
and the Number of Rejected Messages 


Equatorial Array 

Total 

Number 

Messages 

i 

Rejected 

(Elevation) 

Direction Cosine 1 
Rejected Rejected 

(Discrepancy) (S.D.) 

Direction Cosine m 
Rejected Rejected 
(Discrepancy) (S.D.) 

Fort Myers 

52 

20 


1 

1 

Quito 

7 




3 

Lima 

5 





Santiago 

9 


1 


1 

New Foundlancf 

99 

5 

1 


3 

Winkfield 

41 


9 


3 

Johannesburg 

6 





Blossom Point 

73 

15 , 

1 



College 

69 

6 

3 

1 

3 1 

Mojave 

79 

25 




Grand Forks 

59 

4 

1 



Woomera 

6 


1 


1 

Total 

505 

75 

17 

2 

12 3 


Rejection Criteria 

Elevation. All messages the mean elevation of which was below 25°. 

Discrepancy. All direction cosines having a residual error in excess of 

20 Minitrack Counts. 

S.D. All direction cosines having an error with a standard deviation 

in excess of 0.00050, the computed standard deviation being 
based on all measurements within a Minitrack message. 

Rejected measurements are only listed under one heading even if they should have been 
rejected for more than one reason. Thus, the total number of measurements less the 
number of rejected measurements equals the number of measurements used in the 
computations (Table 3) . 
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Table 12 (Continued) 

Total Number of Minitrack Messages 
and the Number of Rejected Messages 



Total 


Direction 

Cosine 1 

Direction 

Cosine m 

Polar Array 

Number 

Messages 

Rejected 

Elevation 

Rejected 

Discrepancy 

Rejected 

(S.D.) 

Rejected 

Discrepancy 

Rejected 

TsTdTT 

Fort Myers 

72 

18 


1 


2 

Quito 

19 


1 

1 

1 


Lima 

21 



2 



Santiago 

22 

1 

4 

1 

3 

2 

New Foundland 

74 

19 

1 




Winkfield 

8 




1 


Johannesburg 

20 


1 




Blossom Point 

82 

24 

1 


2 


College 

76 

40 

1 

6 

1 

6 

Mojave 

70 

10 

1 

1 

1 

1 

Grand Forks 

51 

11 

1 



1 

Woomera 

20 


4 


2 

— 

Total 

535 

123 

15 

12 

11 

12 


Rejection Criteria 

Elevation. All messages the mean elevation of which was below 25°. 

Discrepancy. All direction cosines having a residual error in excess of 

20 Mini track Counts. 

S.D. All direction cosines having an error with a standard deviation 

in excess of 0.00050, the computed standard deviation being 
based on all measurements within a Minitrack message. 

Rejected measurements are only listed under one heading even if they should have been 
rejected for more than one reason. Thus, the total number of measurements less the 
number of re jeered measurements equals the number of measurements used in the 
computations (Table 3). 
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Table 13 

Total Number of Measurements end Number of Rejected Messages 



Total 

Rejected 

Rejected 

Rejected 

Used In 


Measurements 

(Elevation) 

(Discrepancy) 

(S.D.) 

Computation 

Number 

2080 

396 

55 

29 

1600 

Percent 

100 

19.0 

2.6 

1.4 

77.0 


Rejection Criteria 

Elevation. All measurements at a mean elevation below 25°. 

Discrepancy All measurements having a residual error in excess of 20 

Minitrack Counts 

S.D. All measurements having an error with a standard deviation in 

excess of 0.00050, the computed standard deviation being based 
on all measurements within a Minitrack message. 

Rejected measurements are only listed under one heading even if 
they should have been rejected for more than one reason. 
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Table 14 



Residual Measurement 

Errors and Prepass Calibration 

Constants 




for the Woomera Polar Station 






Residual Measurement 






Error 


KC-KSl + KS2 

Arc 

Date 

Time (hrs) 

1 

m 


m 

1 

1.2.66 

5.2 

11 

3 

345 

557 

1 

1.3.66 

18.7 

-10 

5 

344 

545 

1 

1.4.66 

18.8 

- 2 

30 

343 

551 

1 

1.5.66 

18.8 

0 

5 

342 

570 

2 

1.11.66 

17.2 

- 9 

-14 

345 

584 

2 

1.12.66 

17.2 

5 

- 8 

343 

584 

2 

1.13,66 

17.3 

31 

-16 

305 

576 

3 

1.28.66 

14.2 

65 

7 

152 

546 

3 

1.28.66 

22.8 

98 

- 2 

124 

544 

3 

1.29.66 

14.3 

-34 

65 

224 

545 

3 

1.29.66 

22.9 

-16 

12 

222 

541 

3 

1.30.66 

22.9 

3 

4 

221 

540 

3 

2.3.66 

12.6 

14 

10 

160 

541 

4 

3.12.66 

4.8 

1 

-11 

561 

719 

4 

3.12.66 

13.4 

- 4 

0 

561 

717 

4 

3.13.66 

4,9 

0 

- 1 

560 

717 

4 

3.13.66 

13.5 

4 

- 5 

560 

717 

4 

3.15.66 

13.6 

0 

- 1 

560 

717 

4 

3.16.66 

13.7 

- 2 

1 

561 

713 

4 

3.19.66 

3.3 

2 

-12 

563 

716 


100 



Table 15.1 


Comparison of a GEQS-A Arc Determined from Minitrack Data with an Orbit 

Determined from Optical Data 


Start Time 1.1.66 6 hours 

Stop Time 1.5.66 19 hours 

Length of Arc 4.5 days 

Run 1 

Run 2 

Run 3 

RMS Position Difference (m) 

67 

69 

149 

RMS Velocity Difference (m./sec) 

.044 

.046 

0.100 

Maximum Position Difference (m) 

103 

122 

237 

Maximum Velocity Difference (m/sec) 

.063 

.074 

0.140 


Run 1. Minitrack orbit with measurement bias, scale factor, and rotation 
recovery. No ionospheric corrections. 

Run 2. Minitrack orbit with measurement bias and scale factor recovery. 

No ionospheric corrections. 

Run 3. Minitrack orbit with measurement bias, scale factor, and rotation 
recovery. Ionospheric corrections. 

It is of interest to compare the above results with those quoted in "Intercomparison 
of the Minitrack and Optical Tracking Networks using GEOS-I Long Arc Orbital 
Solutions" (by J. G. Marsh, C.E.Doll, R. J. Sandifer and W. A. Taylor, NASA-TND- 
5337, February 1970). For an arc covering almost the identical period, they obtained 
an RMS position difference of 165 meters in a similar comparison. However, they 
did not attempt to recover any station calibration parameters. They made no 
ionospheric co r rections. 
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Table 15.2 


Comparison of a GE QS-A Arc Determined from Minitrack Data with an Orbit 

Determined from Optical Data 


Start Time 1.8.66 0 hours 

Stop Time 1.14.66 5 hours 

Length of Arc 6.2 days 



Run 1 

Run 2 

Run 3 

RMS Position Difference (m) 

86 

116 

160 

RMS Velocity Difference (m/sec) 

.075 

.099 

.127 

Maximum Position Difference (m) 

127 

180 

262 

Maximum Velocity Difference (m./sec) 

.120 

.158 

.180 


Run 1. Minitrack orbit with measurement bias, scale factor, and rotation 
recovery. No ionospheric corrections. 

Run 2. Minitrack orbit with measurement bias and scale factor recovery. 
No ionospheric corrections. 

Run 3. Minitrack orbit with measurement bias, scale factor, and rotation 
recovery. Ionospheric corrections. 
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Table 15.3 


Comparison of a GEQS-A Arc Determined from Minitrack Data with an Orbit- 

Determined from Optical Data 


Start Time 1.28.66 2 hours 

Stop Time 2.4.66 0 hours 

Length of Arc 6.9 days 



Run 1 

Run 2 

Run 3 

RMS Position Difference (m) 

124 

114 

182 

RMS Velocity Difference (m/sec) 

.103 

.099 

.155 

Maximum Position Difference (m) 

194 

181 

272 

Maximum Velocity Difference (m/sec) 

.155 

.153 

.231 


Run 1. Minitrack orbit with measurement bias, scale factor, and rotation 
recovery. No ionospheric corrections. 

Run 2. Minitrack orbit with measurement bias and scale factor recovery. 
No ionospheric corrections. 

Run 3. Minitrack orbit with measurement bias, scale factor, and rotation 
recovery. Ionospheric corrections. 
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Table 15.4 


Comparison of a GEOS-A Arc Determined from Minitrack Data with an Orbit 


Determined from Optica! 

Data 



Start Time 3.12.66 3 hours 

Stop Time 3.19.66 20 hours 

Length of Arc 7. 6 days 

Run 1 

Run 2 

Run 3 

RMS Position Difference (m) 

168 

179 

249 

RMS Velocity Difference (m/sec) 

.142 

.152 

.207 

Maximum Position Difference (m) 

251 

250 

406 

Maximum Velocity Difference (m /sec) 

.207 

.215 

.315 


Run 1. Minitrack orbit with measurement bias, scale factor, and rotation 
recovery. No ionospheric corrections. 

Run 2. Minitrack orbit with measurement bias and scale factor recovery. 
No ionospheric corrections. 

Run 3. Minitrack orbit with measurement bias, scale factor, and rotation 
recovery. Ionospheric corrections. 
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2.6 RECOMMENDED PROCEDURE FOR REDUCING MINITRACK DATA 
2.6.1 Reduction Procedure 


The purpose of this section is to provide a step-by-step procedure to be 
followed for the reduction of Minitrack data and optical data. This procedure 
can be simplified into three or four parts. 


PRE PROCESS 
MINITRACK 
DATA 


PREPROCESS 
OPTICAL 
DATA 


NAP-II 
PROGRAM 
REDUCTION 


REPORT 

RESULTS 


The preprocessing of the Minitrack data and the Optical data can be done 
as a series or parallel effort. The same applies to the NAP-II reduction of the data. 


The steps to a complete reduction of the Minitrack data are: 


1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 


Extract desired data from Minitrack and/or Optical data tapes. 

Sort Minitrack messages in time sequence (not required for 
Optical data when GEOS formatted data tapes are used). 

Preprocess Minitrack and/or Optical data. 


Set up NAP-II control cards and write cards on magnetic 
tape . 


Make PRENAP card updater run. 
NAP-II program execution. 
POSTNAP program execution. 


Iterate on these 
3 steps unti I 
convergence 
achieved 


Report results after reduction has converged. 


A flow diagram for the preprocessing of the Minitrack data is given in 
Figure 1 . Note: In setting up these runs, refer to Appendix A-4 for restriction 
place on card set-up. 
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FIGURE 1 

MINITRACK PREPROCESSING 


Same 

Format 
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The process starts with a magnetic tape containing Minitrack messages from 
several satellites. This tape is input to a program to extract just the messages from 
GEOS-1 satellite and writes them on another magnetic tape in the same format. This 
tape is then input to a program that sorts the messages, and also edits duplicate messages 
and puts them in time sequence on magnetic tape. This tape is now in the right form 
for the preprocessor program. The time-sequenced data tape, along with preprocessor 
control cards, are now input to the Minitrack preprocessor program. The preprocessor 
makes known corrections to the data and converts the corrected phase differences to 
direction cosines. The program outputs a data tape in a format acceptable to NAP-II 
and also a listing of the corrected data . The preprocessor also outputs a bulk of the 
NAP-II control cards dealing with the station times for the data (Category 201, 202, 
and 999). These are used in selection of data to be processed. This ends the pre- 
processing phase of the procedure. The preprocessing of Minitrack data is discussed 
in Section 2.2 of this report. 

Figure 2 is a flow diagram of the preprocessing of optical data. The data input 
to the Optical preprocessing program must be in time sequence. If there are more 
than one tracking systems data to be processed, the data tapes must be "merged" 
onto one data tape. The output from the merge program is a data tape, in GEOS 
format, and a listing of the data. This data tape is input to the Optical preprocessor 
with control cards that specify the stations and times of the data to be processed. 

In the case of the SAO stations, these cards are also used to correct the observation 
times from A . 1 to UTC. The Optical preprocessor outputs a data tape in a format 
acceptable to NAP-II and punched cards used in NAP-II to control data times. 

There is also a printed output of the reformatted and corrected data . The optical 
data preprocessing is discussed in Section 2.3 of this report. 
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FIGURE 2 

OPTICAL PREPROCESSING 
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Once the data has been preprocessed, it is ready to be reduced using the 
NAP-1 1 program. This process is flowed in Figure 3. The first step is to set up the 
NAP-1 1 control cards as specified by the NAP-1! User's Guide. The timing cards 
(Category 201, 202, and 999) are output as a result of the preprocessor. After the 
cards are set up, they are written on magnetic tape, for ease of handling. The 
cards on tape can be changed via the PRENAP card updater. This program is dis- 
cussed in Section 2.4. 1 of this report. 

The NAP-1 1 control cards (now on tape) and the data tape output from the 
preprocessor are now input to the NAP-il program. At the end of each iteration, 
NAP-II outputs a disk file (on tape) with the current values of the error model param 
eters and the measurement discrepancies. Also output is a "single point" data tape 
which can be used on successive iterations to reduce the computation time. When 
using the "single point" data tape for Minitrack data, delete the Category 704 cards 

The output from the NAP-II program is then input to the POSTNAP program 
which, in effect, gives a status report of the reduction and prepares information for 
the next iteration . 

The PRENAP, NAP, POSTNAP sequence is continued until the solution has 
converged. 



NAP- II Con fro I 
Cords 


FIGURE 3 

NAP-II REDUCTION 


Program to 
Write Cards 
on Magnetic 
Tape 


NAP-II 

"Cards" 

Tape 



PRENAP Card 
Updater 

'"NAP-II Update 
Cards 





NAP-II 

Program 


'Updated 
NAP-II 
"Cards" 
V Tape 


Use the Computed 
Single Data Point 
Per Pass 
(FJ37F003) 


FT37F001 

FT37F002 

FT37F003 


POST NAP 
PROGRAM 


Summary 

of 

Results 


/SequencedX 
f Cards With 


\ "Edit" Cardy 
\^Cat.701i 

, 

1 

1 


Data Points / 


Iterations 














2.6.2 Recommended Error Model 

The results of this study indicate the following error models to be sufficient 
for reducing Minitrack data. 

The error model terms to be used are dependent on whether the data was 
collected in the equatorial or polar modes. 

EQUATORIAL ARRAY 
t = lo + b x + r e m 0 + t£ 

m = m 0 + b 3 + s e m 0 + rm 

POLAR ARRAY 
■L — -to + b 3 + s p 4-0 + t4- 

m= m 0 + b 4 + r p £o + rm 

where, 

-to and m 0 are Minitrack measurements in 
Minitrack counts/1000. 

bi, b 2 , b 3 and b 4 are zero-set biases 

r e and r p are rotation terms about the local vertical 
of equatorial -6 and polar m measurements 
respectively 

s g and s p are scale factors of equatorial m and polar 4- 
measurements respectively 

r is the timing error associated with the station 

4, and m are computed rates obtained from the orbital 
equations within NAP-II 

These error model terms are coded for input to NAP-II (see Reference 4, 
The NAP-II User's Guide, Appendix IV-B). 
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The error model terms would correspond to the following term numbers for 
direction cosine data types (measurement code 6 and 7). 

Error Model Term Error Model Term No . 


bi, b 3 , b 3/ b 4 
r e and r p 

s e and Sp 
T 


10 

15 /see changes to MESOLD 
\ in Section 2.4.2 

14 

11 


Note: Consult NAP-1! User's Guide, Appendix IV-2. 
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2.6.3 Typical Set-up for NAP-II Program 

The following is an example of a set of NAP-II control cards. This listing 
is output from the NAP-II program on each run. 

Note that each station is defined twice (category 301, 302, and 303); once 
for the equatorial array and once for the polar array. Following the suivey cards 
are the error model cards (category 601). For ease of recognition, the labels have 
been coded as to which term it is, i.e. FTMYRELB is the FTMYR station equatorial 
"V' measurement zero-set bias "B 11 . FTMYRPLB is the polar "-C" measurement bias. 
Other codes used are "ELS" and "PLS" for equatorial "l" scale and polar scale, 
"ELR" and "PLR" for equatorial rotation and polar rotation. 

It is recommended that a new user of NAP-II take the following listing and, 
using the NAP-II User's Guide, cross-reference the codes used in the listing to 
their meaning given in the User's Guide. 

See Appendix A-4 for restriction on deck setup. 
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APPENDIX A-l 

CONSTANTS USED IN NAP-II 


A— 1-1 



A- 1.1 


THE GRAVITY MODEL USED 


The gravity model used for data processing under this contract v/as the Smithsonian Ml 
model, (Reference 1), which is built into the NAP program, modified by the GEOS-A 
resonant harmonics (Reference 2). The GEOS-A resonant harmonics were not obtained 
from the original source, but from (Reference 3). 

The gravitational mass of the Earth was taken as .3986032 D15 rn 3 /sec 2 , which is the 
NAP default value . 

The Smithsonian Ml values for the spherical harmonics coefficients are listed below: 


(N, M) 

C(N,M) 


S(N,M) 


(Q, 0) 

1. 


0. 


Or 0) 

0. 


0. 


(2, 0) 

- .108264500002 

D-2 

0. 


(3, 0) 

.254599999999 

D-5 

0. . 


(4, 0) 

. 164900000002 

D-5 

0. 


(5, 0) 

.21 

D-6 

0. 


(6, 0) 

- .645999 999993 

D~6 

0. 


(7, 0) 

.332999999998 

D-6 

0. 


(8, 0) 

' .269999999997 

D-6 

0. 


(9, 0) 

.529999999998 

D-7 

0. 


(10, 0) 

.540000000002 

D-7 

0. 


(11, 0) 

- .302 

D-6 

0. 


(12, 0) 

.357000000006 

D-6 

0. 


(13, 0) 

.113999999999 

D-6 

0. 


04, 0) 

- .178999999999 

D-6 

0. 


Or 1) 

0. 


0. 


(2/D 

0. 


0. 


(3, 1) 

.209111899862 

D-5 

.287312837632 

D-6 

(4, 1) 

- .542646846485 

D-6 

- .444932466779 

D-6 

(5r 1) 

- .676515582477 

D-7 

- .88203930748 

D-7 

<6r 1) 

- .36979402246 

D-7 

- .21243486397 

D-7 

(7r 1) 

.14418923578 

D-6 

.114180308536 

D-6 

(8r 1) 

- .515388203208 

D-7 

.446669776114 

D-7 

(9, 1) 

.76024995806 

D-7 

.779743547594 

D-8 

(10, 1) 

.648810099708 

D-7 

- .778572119639 

D-7 

m -\\ 
\ • • / • / 

m9Q70fin7onn 

• v/ Mk.'/ £.w/ / >jvv 

rV: "7 
t-v / 

OOC/OOA70A 4 A 
. <J/V/ 

o 

u~ o 

02, 1) 

- .922805782847 

D-7 

- .401958347128 

D-7 

(13, 1) 

0. 


0. 


04, 1) 

- .788307409236 

D-8 

.278535284599 

D-8 


A-l-2 


(N,M) 

C(N,M) 


S(N,M) 


(2, 2) 

.153563789676 

D-5 

~ .872066750118 

D-6 

(3, 2) 

.250703728736 

D-6 

- .183761983737 

D-6 

(4, 2) 

.737902432577 

D-7 

.147804093315 

D-6 

(5, 2) 

.102117707526 

D-6 

- .375456547484 

D-7 

(6, 2) 

.858383031385 

D-8 

- .455316216648 

D-7 

(7, 2) 

.362552678174 

D-7 

.162351886106 

D-7 

(8, 2) 

.21354899813 

D-8 

.320323497197 

D-8 

(9, 2) 

- .27706980487 

D-9 

.242436079257 

D-8 

(10, 2) 

- .624317809516 

D-8 

- .249727123811 

D-8 

(11,2) 

0 . 


0 . 


(12, 2) 

- .469893640425 

D-8 

-• .232665783121 

D-9 

(3, 3) 

.782277124813 

D-7 ' 

.225898207164 

D-6 

(4, 3) 

.508569773273 

D-7 

- .113546717887 

D-7 

(5, 3) 

- .171778742838 

D-7 

.231240615425 

D-9 

(6, 3) 

- .111963004094 

D-8 

.642750579065 

D-9 

(7, 3) 

.352147606138 

D-8 

.253546276422 

D-9 

(8, 3) 

- .374070918995 

D-9 

.404400993503 

D-10 

(9, 3) 

0 . 


0 . 


(10, 3) 

- .378977225408 

D-9 

.174912565571 

D-9 

(4, 4) 

- .111982938 75 

D-8 

.485963696475 

D-8 

(5, 4) 

- .206336328204 

D-8 

.498321698304 

D-9 

(6, 4) 

- .166560812946 

D-9 

- .196087502512 

D-8 

(7, 4) 

- .322776554575 

D-9 

- .216600056362 

D-9 

(8, 4) 

- .276702101966 

D-9 

- .156623831302 

D-10 

(9, 4) 

0 . 

• 

0 . 


(10, 4) 

- .435831320535 

D-10 

- .653746980798 

t>-10 

(5, 5) 

.384103946188 

D-9 

- .145764420605 

0-8 

(6, 5) 

- .252611520257 

D-9 

- .369636026451 

D-9 

(7/ 5) 

.268980462153 

D-10 

.191117696793 

D-10 

(8, 5) 

• - .959291939079 

D-n 

.213578205308 

D-10 

(6, 6) 

- .931919036559 

D-l 1 

- .36111862666 

D-10 

(7, 6) 

- .145066373826 

D-10 

.437281413924 

0-11 

(8, 6) 

- .474787418211 

D- 12 

.88813175875 

0-11 
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S(.N>M) 


(N,M ) C(N,M) 


(7, 7 ) 

.102027945339 

on 

.178085140954' 

on 

(8, 7) 

- .443617740607 

D- 13 

.158070689181 

012 

(8, 8) 

1 

- .316141378363 

D- 12 

.130025889487 

012 

(12, 12) 

0 . 


0 . 


(13, 12) 

- .1082186306 

018 

.895602460152 

019 

(14, 12) 

0 . 


0 . 


(15, 12) 

- .114596608272 

019 

. .107203278705 

019 

(13, 13) 

- .274440770475 

D- 19 

.365921027295 

020 

(14, 13) 

0 . 


0 . 


(15, 13) 

- . 127051846403 

020 

- .133101934329 

020 

(14, 14) 

0 . 


0 . 


(15, 14) 

.219788278315 

022 

- .532258360739 

022 


The GEOS-A resonant spherical harmonics coefficients are listed below: 


(N,M) 

C(N,M) 


S(N,M) 


(13,12) 

-.126299 

018 

.165220 

018 

(14, 12) 

.139978 

020 

-.131772 

0-19 

(15,12) 

-.138126 

019 

-.189639 

020 
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A-!. 2 THE EARTH MODEL USED AND STATION COORDINATES 

The C-5 Earth Model was adopted for the data processing. According to this model 
the Earth is defined geometrically by the following constants: 


Rotation Rate 
Semi -Major Axis 
Eccentricity Squared 

= .7292115854937 D-4 radians/sec 

= 6378165 meters 

= .669454185459 D-2. 


The station 

coordinates are 

given below: 





Station Coordinates (Optical) 


Station 

Station ID 

Latitude 

Longitude 

Height Above Geoid 

COLDLK 

29424 

54°44 , 37.26" 

249°57'21 .90" 

548 meters 

NEWFL 

31032 

47°44‘28.73" 

307°16 , 46.67" 

58 meters 

JUPTH 

37073 

27° 1 '14.33" 

279°53 , 12.72" 

- 41 meters 

JUPTR 

29010 

27° 1 '14.23" 

279°53'12.95" 

36 meters 

JBC4L 

37074 

27° 1 '14.55" 

279°53'12.76" 

- 38 meters 

COLEG 

31033 

64°52'17.78 " 

212°9 I 37.29" 

139 meters 

OOMER 

31024 

-31°23 i 26.96" 

136°52'14.25 " 

148 meters 

EDINB 

37036 

26°22'46.35 " 

261°40'7.34" 

15 meters 

JUP40 

37072 

27° 1 ’14.39" 

279°53'12.49 " 

38 meters 

GFORK 

31034 

48° 1 '20.81 " 

262°59 , 19.55" 

200 meters 

ROSMA 

31042 

35°12'7.03" 

277°7'40.81 " 

857 meters 

PURIO 

37040 

18°15'28 .30" 

294°0'23.63" 

5 meters 

ORGAN 

29001 

32°25'24.70 " 

253 0 26'48.29" 

1610 meters 

GSFCP 

37043 

39° 1 '14.78" 

283°10'20.39" 

1 meter 

BPOIN 

31021 

38°25'49.44" 

282°54 , 48.65" 

- 50 meters 

EDWAFB 

29425 

34°57'50. 17" 

242°5'7.80" 

754 meters 
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Station Coordinates (Optical) (Continued) 


Station 

Station ID 

Latitude 

Longitude 

Height Above 

JUP24 

37071 

'27° 1' 14. 00" 

279°53'12.30" 

- 38 meters 

DENVR 

3 7045 

39°38 , 47.54" 

255 0 23'38.52" 

1751 meters 

COLBA 

37037 

38°53'35 .81 " 

267°47 , 40.85" 

218 meters 

FTMYR 

31022 

26°32'53.08" 

278°8'3.80" 

- 42 meters 

BERMD 

37039 

32°21 ’48.94" 

295°20'34. 18" 

- 28 meters 

MOJAV 

31030 

35°19 , 47.57" 

243°5 '59. 18" 

874 meters 

SATAG 

31028 

-33°8'58 .76" 

289° 19 '52. 59" 

705 meters 
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Station Coordinates (Minitrack) 


Station 

Latitude 

Longitude 

Height Above Geoid 

Fort Myers 

26° 32 '5 3. 08" 

278°8 , 3.80" 

- 42 meters 

Quito 

-37 ‘22. 63" 

231°25'15 .23" 

3554 meters 

Lima 

-1 1°46'37.56" 

282°50'58.86" 

34 meters 

Santiago 

-33°8'58.76" 

289°19‘52 .59” 

705 meters 

New Foundland 

47°44'28.73" 

307° 16 '46. 67" 

58 meters 

Winkfield 

51°26'40 .67" 

359° 18 ‘8. 35" 

76 meters 

Johannesburg 

-25°53'2.7" 

27° 42 ‘25 . 4 1 " 

1546 meters 

Blossom Point 

38° 25 '49. 44" 

282°54'48.65" 

- 50 meters 

College 

64°52 , 17.78" 

212 0 9'37.29" 

139 meters 

Mojave 

35°19'47.57" 

243°5 '59.18 " 

874 meters 

Grand Fork 

48°1 *20.81 " 

262 0 59'19.55" 

200 meters 

Woomera 

-31°23 , 26.96" 

136°52' 14.25 " 

148 meters 
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APPENDIX A -2 

SPECIAL PROGRAMS WRITTEN TO AID 
MiNITRACK DATA REDUCTION 
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A -2. 1 PROGRAM FOR EXTRACTING MINITRACK MESSAGES 
FOR GEOS 1 AND SAMPLE JCL 



PROGRAM FOR EXTRACTING M NI TRACK MESSAGES FOR SATELLITE ID 65891 
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A PROGRAM FOR MERGING TWO OPTICAL DATA TAPES 
(IN A GEOS FORMAT) INTO A SINGLE TAPE - SAMPLE JCL 
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Appendix A-4 

Constraints Imposed on 
NAP-II Deck Setup 
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A-4.1 


CONSTRAINTS IMPOSED BY MINITRACK PREPROCESSOR 


There are constraints imposed on the NAP-II deck setup by the preprocessing pro- 
gram . This appendix is an attempt to outl ine the const rai nts in a concise , quick- 
reference form . 

The Minitrack Preprocessor (Section 2.2.3) imposes an order to the stations in NAP-II 
(Key 1, Category 201), depending on the order the Minitrack data is preprocessed 
(data set 5 of Minitrack Preprocessor). The station ID's required by the NAP-11 
program are generated from the fifty-fifth and fifty-sixth characters of the Minitrack 
message format (Section 2.2.1). This code is used for the equatorial data, and the 
code plus 100 is used for the polar data. 

Example: Station constants for Minitrack processor are in the following order: 

Ft. Meyers 

Quito 

Lima 

Santiago 

N ewfoundland . 

The following information would be generated by the preprocessor for each station: 

Equatorial Polar 



Sta . ID 

Sta . No . 

Sta. ID 

Sta . No 

Ft . Meyers 

03 

1 

103 

2 

Quito 

05 

3 

105 

4 

Lima 

06 

5 

106 

6 

Santiago 

08 

7 

108 

8 

Newfoundland 

12 

9 

112 

10 
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However,, if the station's constants have been in order of: 
Newfoundland 
Santiago 
Lima 
Quito 

Ft. Meyers, 

then the information generated would be 



Equatorial 


Polar 


Sta . ID Sta . No . 

Sta. ID 

Sta . No 

Newfoundland 

12 

1 

112 

2 

Santiago 

08 

3 

108 

4 

Li ma 

06 

5 

106 

6 

Quito 

05 

7 

105 

8 

Ft. Meyers 

03 

9 

103 

10 


The station numbers are required on the NAP-II cards (Key 1 of Category 301 and 
Key 2 of Category 201 and 202) punched by the preprocessor program. The station 
identification is required in Keys 3 and 4 of the Category 301 card. 

Pass numbers are also generated in the preprocessor and are punched in Key 3 of the 
Category 201, 202, and 999 cards. 

In summary, the following information is generated for the Minitrack Preprocessor for 
the NAP-11 program. 


A. STATION ID 

Category 30 1 

Key 4 

B. STATION NUMBER 

Category 301 

Key 1 

1) 

Even number for polar 



2) 

Odd number for 
equatorial 
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C. Start and Stop times of the 

data pass Category 201 and 202 

D. End of pass control card Category 999. 
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CONSTRAINTS IMPOSED BY THE PRE-NAP PROGRAM 


This program is used to modify, add, or delete NAP-1! control cards. 

The following are the instructions for its usage: 

INS 1 1 Inserts cards after sequence number 1 1 . 

DEL II 12 Deletes cards between (and including) sequence numbers II and 12. 

REP II 12 Deletes cards between (and including) sequence numbers 11 and 

12, and then inserts cards 

MOD II 12 Modifies cards between (and including) sequence numbers II 
and 12 . 

Format for above four cards (A4, x, 15, 5x, 15) note that INS, DEL, REP, and MOD 
start in column 1 and that column 2 should be blank. 

INS and REP must be followed by the cards to be inserted. The cards must be fol- 
lowed by a blank card to signify the end of the particular set of inserted cards. 

The MOD cards must be followed by a modifying card. Columns not to be modified 
should be left blank. 

When using the computed mid-point Minitrack data, the Category 704 cards have to 
be deleted from the control cards. The "704" cards should be used only when the 
entire data spans are used. 
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A-4.3 


CONSTRAINTS IMPOSED ON NAP-1 1 SETUP 


A very basic assumption made when processing Minitrack data is that a good estimate 
of the state vector is available. This is the basis ->f the whole process. 

The a priori sigma information should be realistically assigned with respect to the 
error model terms and the observed data . 

When processing with the "single-data, point pen pass" data file, be sure to delete 
Category 704 cards that were used for the processing of the entire observed data 
file . 

Measurement numbers (Category 701) have a limitation as to their assignment. Odd 
measurement numbers are reserved for cosine alpha (£) measurements, and even 
numbers are reserved for cosine beta (m) measurements. This limitation is imposed 
by the Post-NAP processor. 

Error model term definition cards (Category 601 cards) must be in order of totally 
stable, arc stable, and pass stable. There can be only one 601 card for each error 
model term. Duplicates inserted to override a previous 601 card are not permitted. 
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ADDITIONAL CONSTRAINTS 


This section is reserved for constraints not mentioned above that have not come to 
light. If, in the course of reducing Minitrack dura, using the procedure outlined 
in this report, new constraints are found, please note them in this section. 
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